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SUMMARY 


This  final  report  on  Contract  F30602 -80-C-0042  summar¬ 
izes  the  AIL  investigation  to  determine  and  demonstrate  suitable 
algorithms  and  corresponding  computer  software  in  order  to  pro¬ 
vide  open-loop  adaptive  control  of  the  Microwave  Transversal 
Equalizer  (MTE)  previously  developed  by  AIL  under  Contract  F3062- 
78-C-0352. 

The  program  task  objectives  have  been  successfully  accomp¬ 
lished  by  AIL  with  the  application  of  Fast  Fourier  Transform  (FFT) 
techniques  to  provide  an  algorithm  programmed  on  the  RADC  HP  2100A 
computer.  This  procedure  will  ultimately  be  utilized  to  determine 
the  MTE  amplitude  and  time  delay  adjustments  necessary  to  equalize 
the  distortion  introduced  by  an  arbitrary  network  in  series  with 
the  equalizer. 

Verification  of  the  developed  FFT  and  associated  software 
program  was  accomplished  at  AIL  with  a  DEC-20  computer  using 
artificially  generated  simulated  time  sidelobe  distortion  typical 
of  the  expected  range  of  levels.  Translation  of  the  algorithm 
from  the  DEC-20  computer  format  to  the  HP  2100A  computer  format 
has  been  completed,  and  verification  of  the  software  program  using 
RADC  I  and  Q  data  is  underway.  Preliminary  data  suggests  that  the 
program  translation  task  has  been  successfully  accomplished  by  AIL. 
It  is  expected  that  optimal  data  acquisition  techniques  and/or 
preferred  computer  interface  procedures  will  be  defined  as  a 
result  of  related  post-delivery  activities  at  RADC. 
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1.  INTRODUCTION 


The  objective  of  Contract  F30602-80-C-0042  is  to  develop 
an  algorithm,  and  relevant  computer  software  programming  for  the 
RADC  HP  2100A  computer,  to  facilitate  open-loop  control  of  the 
MTE  in  order  to  minimize  time  sidelobe  distortion.  The  result¬ 
ant  algorithm  using  FFT  techniques  and  associated  programming 
has  been  validated  with  artificially  generated  simulated  distor¬ 
tion  data  using  the  AIL  DEC-20  computer  and  the  HP  2100A  computer 
at  RADC.  A  similar  post-delivery  validation  effort  will  be  con¬ 
ducted  by  RADC  personnel  using  real-time  distortion  data  and  tne 
HP  2I00A  computer. 

Section  2.0  will  contain  a  brief  discussion  of  MTE  opera¬ 
tional  fundamentals  in  order  to  establish  an  inter-relationship 
with  the  algorithm  development  program.  Section  3.0  will  pre¬ 
sent  a  detailed  treatment  of  the  FFT  techniques  applied  to 
algorithm  development  for  the  MTE.  Section  4.0  will  discuss 
the  resultant  computer  program  including  applicable  operating 
and  test  procedures.  Finally,  Section  5.0  will  present  con¬ 
clusions  and  recommendations  for  a  logical  continuation  of 
effort  required  to  provide  closed-loop  operation  of  the  MTE. 
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2.  MTE  TECHNICAL  DISCUSSION 


A  discussion  of  distortion  considerations  and  the  funda 
mentals  of  the  MTE  operation  are  presented  in  order  to  provide 
the  technical  basis  for  computer  control  and  adaptive  equaliza¬ 
tion  of  resultant  time  sidelobe  distortion.  A  detailed  descrip¬ 
tion  of  the  MTE,  including  operating  procedures,  is  included  in 
Final  Report  RADC-TR-80-1 21  of  31  January  1980  titled  "Microwave 
Transversal  Equalizer". 


2.1  PAIRED  ECHO  THEORY  AND  DISTORTION  CONSIDERATIONS 

Algorithms  for  determining  the  settings  of  attenuation 
and  phase  for  each  loop  of  a  transversal  equalizer  rely  heavily 
on  the  paired  echo  concept.  Accordingly,  a  brief  review  of  this 
theory  is  given  below  stressing  its  applicability  to  the  present 
program . 


2.1.1  Paired  Echo  Theory 

The  paired  echo  concept  and  its  application  to  the 
design  of  a  transveral  equalizer  (MTE)  was  described  in  Ref.  1 
It  should  be  noted  that  this  work  wis  supported  by  RADC.  This 
paper  also  gives  a  bibliography  of  prior  work  on  non-microwave 
MTE’s  and  the  now  classic  reference  on  Chirp  radar  by  Klander, 
et  al .  The  latter  emphasized  the  use  of  paired  echo  theory  to 
predict  time  sidelobe  levels  in  Chirp  radars.  The  essential 
aspects  of  the  paired  echo  theory  are  summarized  below. 


Ref.  1  -  J.J.  Taub  and  G.P.  Kurpis,  "Microwave  Transversal  Equal 
izer".  Microwave  Journal,  1969. 
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Signals  which  can  be  characterized  by  a  time  function 
of  its  Fourier  transform  (amplitude  and  phase  vs  frequency  func¬ 
tions),  are  subject  to  distortion  when  propagating  through  a  micro- 
wave  component  or  a  system  of  microwave  components  (such  as  a  Chirp 
radar).  This  distortion  occurs  because  the  system's  transfer  func¬ 
tion  possesses  neither  perfectly  constant  gain  nor  perfectly  linear 
phase  over  the  spectrum  of  the  signal. 

We  typically  represent  the  transfer  function  in  the  fre¬ 
quency  domain  and  need  to  predict  its  effect  on  time  domain  dis¬ 
tortion.  By  using  the  paired  theory  we  can  make  relatively  quick 
conversions  between  the  frequency  domain  and  the  time  domain  and 
vice  versa. 

The  system  transfer  function,  or  frequency  response,  of 
an  aribtrary  system  can  be  defined  as: 

H(a))  *  A(u))  (2-1) 

where  A  and  B  are  respective  gain  and  phase  functions  defined  in 
Reference  1 . 

We  must  now  assume  that  the  signal  transmitted  through 
the  system  has  spectral  components  that  are  band  limited.  This 
is  a  safe  assumption  for  most  systems.  For  example,  in  many  radar 
applications  the  signal's  spectrum  is  confined  to  the  electronic 
bandwidth  of  the  final  power  amplifier.  Within  a  band  limit  of 
u)^  to  we  can  rigorously  represent  the  A  and  B  functions  as 
Fourier  series: 
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In  a  distortionless  network  A(w)  =  and  B(^)  = 

Hence  t'^e  remaining  terms  contribute  to  distortion.  This 
representation  leads  to  the  time  domain  response  including 
distortion  and  so  that  for  a  band  limited  signal 
V.(;.j)  and  its  time  function  equivalent  V^(t)  we  obtain 

V^(t)  .  »,  (T  +  !>  )  +  a  Z  lt*l,  -  ne) 

n  =  I 

„  (^-4) 

*  S  „S,  E„-  »,  <“'>0 

where  ( t )  is  the  output  time  function  and  the  constants  E  ^  and 
0  n  + 

Ep_  represent  echo  or  time  sidelobe  amplitude  levels.  These  l^v/els 

are  related  to  the  Fourier  distortion  coefficients  in  equations 

(2-2)  and  (2-3).  For  the  case  where  0  and  <p  =0  they  reduce  to 

0  0 
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(2-5) 


A  typical  example  of  a  distorted  microwave  pulse  is  given 
in  Figure  2-1.  The  echoes  are  clearly  displayed.  This  analysis 
desmonstrates  that  a  frequency  domain  transfer  function  can  be 
represented  by  a  Fourier  series  which  enables  rapid  calculation 
of  the  time  domain  response.  Conversely,  a  measurement  of  the 
time  domain  response  which  yields  sidelobe  levels  can  be  used  to 
rapidly  calculate  the  a  and  b  coefficients  thereby  yielding  the 
frequency  domain  transfer  function.  The  MTE  settings  cancel 
distortion  by  injecting  equal  and  opposite  echoes  in  cascade 
with  the  system. 

2.1.2  Use  of  Paired  Echo  Theory 

Since  the  MTE  design  and  adjustment  procedures  are  based 
on  cancelling  distortion  echoes,  algorithms  can  be  developed  by 
taking  either  frequency  domain  (I  and  Q  data  )  measurements  or 
time  domain  responses  (time  sidelobe  levels  in  d3)  and  convert¬ 
ing  them  into  the  necessary  MTE  loop  attenuation  and  phase  set¬ 
tings  to  produce  cancelling  echoes.  The  key  point  of  this  dis¬ 
cussion  is  that  use  of  paired  echo  theory  simplifies  computation 
and  therefore  significantly  reduces  computer  time.  Furthermore 
it  provides  the  flexibility  to  work  with  either  frequency  or 
time  domain  data.  In  this  program  the  actual  procedure  used 
was  to  use  I  and  Q  data  (frequency  domain  measurements). 
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Figure  2  1.  Time  Sideiobe  Distortion 


2.2  MICROWAVE  TRANSVERAL  EQUALIZER  OPERATION 

The  MTE  has  been  designed  by  AIL  to  cancel  the  para¬ 
sitic  time  sidelobes  caused  by  distortion,  by  generating  its 
own  artificial  equalizer  sidelobes,  whose  envelopes  are  prop¬ 
erly  delayed  (or  advanced)  to  coincide  in  time  with  the  para¬ 
sitic  sidelobes.  In  addition,  the  amplitude  and  phase  of  each 
artificial  sidelobe  can  be  separately  adjusted  to  produce  an 
equal  amplitude  phase  reversed  replica  of  the  corresponding 
parasitic  sidelobe. 

Generation  of  the  equalizer  echoes  is  accomplished 
by  tapping  off  portions  of  energy  from  the  main  signal  lobe 
on  the  main  transmission  path,  directing  these  portions  through 
properly  adjusted  delay  lines  and  reintroducing  the  delayed 
signal  portions  back  into  the  main  transmission  path  with 
proper  amplitudes  and  polarity.  The  time  domain  output 
response  of  an  ideal  MTE  is  shown  in  Figure  2-2.  It  should 
be  noted  however  that  the  actual  response  will  consist  of  a 
train  of  echoes  with  decreasing  amplitudes  (0.75  dB/tap)  due 
to  the  insertion  loss  of  the  mainline  tapped  delay  elements. 

A  correction  for  the  aforementioned  MTE  echo  transfer  char¬ 
acteristics  will  be  necessary  in  order  to  achieve  the  desired 
tap  attenuator  setting  with  adaptive  closed  loop  operation 
(computer  controlled). 
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2-2.  Time  Domain  Output  of  Microwave  Transversal  Equalizer. 


The  MTE  consists  of  four  cascade  connected  microwave 


integrated  circuit  mainlines  and  32  secondary  lines  for  the 
purpose  of  generating  a  train  of  8  pre-  and  post-mai nl obe  echoes 
(16  total).  A  3.3  nsec  time  delay  between  echoes  of  each  train 
is  provided.*  The  amplitude  level  and  time  delay  of  the  selected 
16  of  the  32  possible  echoes  are  individually  adjusted  with  the 
PIN  diode  variable  attenuator  (electronically  controlled),  a 
line  stretcher  variable  delay  (mechanically  controlled),  included 
in  each  secondary  line.  The  remaining  16  echo  taps  are  terminated, 
however  these  taps  are  available  for  use  as  required. 

It  should  be  noted  that  the  transfer  characteristics  of 
each  of  the  16  variable  attenuators,  as  well  as  the  16  variable 
solid  state  time  delay  units  (when  they  are  developed)  are  required 
to  properly  achieve  adaptive  MTE  control. 


3.  ALGORITHM  DEVELOPMENT 

An  algorithm  has  been  developed  using  FFT  techniques 
to  provide  open  loop  control  of  the  MTE.  The  technique  will 
be  described  in  terms  of  the  relevant  input  and  output  formats, 
as  well  as  the  FFT  algorithm. 

3.1  INPUT  FORMAT 

The  input  to  the  FFT  consists  of  a-set  of  complex  correl¬ 
ations  as  a  function  of  frequency.  These  samples  can  be  expressed 
by 

H(ju)i)  =  I  (o),.)  +  jQ(a3.)  (3-1) 

where 

oj^  is  the  i^  frequency 
I(u)^)  is  the  in-phase  frequency  response 
Q(oj^)  is  the  quad-phase  frequency  response 

For  the  problem  at  hand,  the  following  frequencies  are  considered: 

=  2ff  [f^^  +  (i-1)  Af],  i  =  1  ,2  , .  256  (  3-2) 

where 

fj_  =  3.1  GHz 

Af  =  1 .17647MHz  =  1 .18  MHz 

The  last  frequency  is  found  by  setting  i  =  256;  that  is, 
f^  =  3-1  +  255  X  1.18  X  10’^  =  3.4  GHz 

As  a  result,  a  total  of  256  pairs  of  measurement  data  are  involved 
for  a  bandwidth  of  300  MHz  as  shown  in  the  format  presented  in 
Tabl e  3-1 . 
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TABLE  3-1.  TAPE  FORMAT  FOR  THE  FREQUENCY  RESPONSE 


Word  Index 

Notation 

Frequency 

1 

) 

3.1  GHz 

2 

f 

0(  jo)^  ) 

3.1  GHz 

3 

I  ( 3^2 ) 

3.1  +  Af 

4 

CM 

3 

O' 

3.1  +  Af 

0 

0 

0 

0 

0 

0 

0 

0 

0 

509 

I(j»255) 

3.1  +254  At 

510 

0(3^255) 

3.1  +  254  Af 

51  1 

25^ 

3.4  GHz 

512 

3.4  GHz 

Note : 

1  . 


-  3.4  -  3.1 
'  255  X  10-^ 


1.18  MHz 


2.  Each  word  is  represented  by  8  bits  in  binary  2's  complemert 
notation  ;  i . e . : 

-128  <  1,0  <  +  127 


The  above  representation  defines  the  formula  transforma¬ 
tion  required  from  the  input  medium  to  the  computer. 
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The  input  samples  are  represented  by 
an  8  bit  binary  two's  complement  notation.  Thus,  a  word  is 
defined  to  be  8  bit  with  the  value  of  the  word  ranged  from  127 
to  -128.  In  floating  point  notation  the  value  is  ranged  from 

127/128  =  0.992188  to 
-128/128  =  -1 

3.2  OUTPUT  FORMAT 

The  output  of  the  FFT  consists  of  a  set  of  =  17 
complex  MTE  tap  coefficients  which  represents  a  truncated  Fourier 
series  representation  of  the  equalizer  transfer  function.  The 
total  number  of  relative  equalizer  taps  is  sixteen  with  the  center 
tap  (number  9)  taken  as  the  reference  tap. 

The  output  tap  coefficient  will  be  presented  by  decimal 
representation  with  the  absolute  magnitude  of  the  tap  less  than 
unity,  that  is;  -1  1  <  1  for  i  =  1,2,. ..17  and  phase  angle  is 

expressed  in  degree  units. 

Further  data  format  transformation  will  be  needed  to 
match  the  exact  setting  on  the  MTE  attenuator  dial. 

3.3  FFT  ALGORITHM 

Consider  again  the  measured  frequency  response  of  a  net¬ 
work  with  a  transfer  function  The  measurement  can  be  func¬ 

tionally  described  by  the  block  diagram  shown  in  Figure  3-1.  The 
input  signal  to  the  device  can  be  written  as 

X.(t)  =  A,-cos  (u).t  +  0^)  (3-3) 
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I(jw)  Q(j«) 


FIGURE  3-1.  MEASUREMENT  SETUP 
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where 


.  is  the  in-band  frequency 
0^  is  the  input  phase  angle 
A  is  the  amplitude  of  the  signal 

The  response  of  the  network  is 

Y^(t)  =  cos  (oj^t  +  0.)  (3-4) 

where 

0^  is  the  corresponding  phase  angle  for  the  i^*^  signal 
is  the  corresponding  amplitude  for  the  i^*^  signal 

The  quadrature  hybrid  is  used  to  generate  the  in-phase  and  the 
quad-phase  components: 


Kju:.) 

1 

T 

A.B. 

cos 

cos  (a)^t  +  0^. ) 

dt 

(3-5) 

and 

Q(j..o^) 

1 

T 

0 

T 

/  A.B  . 

si  n 

cos  (a)^.t  +  0^) 

dt 

(3-6) 

0 


where  T  is  defined  as  the  correlation  time, 
can  be  reduced  to 

AB. 

K cos 
AB. 

Q(joj^)  =  —  sin  {0o'®i^ 


The  correlator  output 

(3-7) 

(3-8) 


The  gain  and  phase  response  as  a  function  of  frequency  can  be 
expressed  by 


G,  =  /l(  *  *  I  ■ 

1  q. 

0.  =  tan"'  (j-!-)  +  0  (3-10) 

1  1^0 
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where  the  short  hand  notation  is  adapted;  that  is 
I,-  =  )  and  Q^.  =  0(ja)^. ) 

Since  A  and  0^  are  not  a  function  of  frequency,  these  measurement 
data  represent  the  transfer  function  The  impulse  response 

h(nT)  can  be  found  by  discrete  Fourier  transform  method. 

A  few  assumptions  will  be  made  such  that  the  Fourier  series 
approach  can  be  applied: 

(1)  Signals  passing  through  the  network  will  be  bandlimited  to 
300  MHz.  To  provide  equalization  over  this  bandwidth  requires 

a  corresponding  time  delay.  This  defines  the  delay  between  taps 
to  be: 

T  *  ®  ^ - r  =  3.33  nsec 

300  X  10 

(2)  Since  any  digital  filter  spectrum  to  be  realized  will  be 
made  periodic  of  period  f^,  then  the  spectrum  can  be 
represented  by  a  Fourier  series: 

00 

H{ja))  *  E  Ck  exp  (jkwT)  (3-11) 

k»  -a. 

where  Cj^  is  the  k^^  Fourier  coefficient  defined  by 

c.  ^  ^  H(ja))  exp  (-jku)T)  du)  (3-12) 

■ 

Generally,  if  the  gain  response  is  an  even  function  and  the  phase 
response  is  an  odd  function,  that  is 

G(j(ij)  *  G(-ja))  (3-13) 

0(ju))  *  -0(-joo)  (3-14) 
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then  the  resultant  coefficients  will  be  real.  Note  that  these 
requirements  that  G(j(jj)  and  be  even  and  odd  respectively, 

lead  to  similar  requirements  that  the  in-phase  and  the  quadrature-phase 
components  to  be  even  and  odd,  respectively.  In  practice,  however, 
these  symmetries  cannot  be  assumed  as  shown  in  Figure  3-2.  As  a 
result,  the  coefficients  will  be  comp! ex;  that  is: 


c,^  =  a^  +  j  b^,  k  =  1,2, . ,17  (3-15) 

where 

9k  ■  /  ‘I  *  ol 

9|^  =  tan"’  (b|^/a|^)  (3-17) 


To  achieve  a  set  of  real  equalizer  coefficients  usually 
means  to  double  the  upper  frequency  bound  of  the  signal  samples  (w^). 
When  the  upper  frequency  is  doubled,  the  second  half  of  the  spectrum 
is  copied  from  the  original  bandwidth  with  the  proper  polarity 
attached.  Examples  to  illustrate  these  properties  will  be  presented 
in  the  next  section. 

3.4  FFT  ILLUSTRATIVE  EXAMPLES 

A  few  examples  illustrates  the  Fourier  transform  approach 
to  well  defined  H(ja))  will  be  given  in  this  section. 

3.4.1  Lowpass  Filter 

Figure  3-3  shows  a  specification  for  a  lowpass  transfer 
function  with  cutoff  frequency  set  at  0.5  rad/sec  and  the  upper  fre¬ 
quency  bound  (u)^)  set  at  2  rad/sec.  The  spectrum  between  w^/2  to 
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FIGURE 


G(j6>) 


(a)  GAIN  SPECTRUM 

0(jcj) 


3-2.  GAIN  AND  PHASE  RESPONSE  (GENERAL  CASE  -  ASYMMETRIC) 
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oj^  is  duplicated  to  provide  an  even  function  for  the  in-phase 
response  I(jiii).  The  quad-phase  response  is  set  to  zero. 

The  Fourier  coefficients  can  be  readily  found  to  be 

c,^  =  sin  (k^),  k  =  -N,  -N  +  1,...N  ^2 

^  k  0  and  N  =  8 


and  Cq  =  0.5. 

For  the  case  N  =  8,  or  17-tap  equalizer  the  response  of 
the  network  can  be  plotted  as  shown  in  Figure  3-4. 

Generally,  the  performance  of  such  a  filter  indicates 
a  passband  ripple  of  0.75  dB  and  a  stopband  rejection  of  approxi¬ 
mately  -20  dB.  Direct  truncation  of  the  Fourier  series  leads  to 
the  well-known  Gibbs  phenomenon.  The  ripple,  in  generaly,  cannot 
be  reduced  by  simply  including  more  taps  as  displayed  by  the  25 
tap  version  shown  in  Figure  3-5. 

The  reduction  of  the  passband  ripple  and  stopband  attenu¬ 
ation  is  later  approached  by  finding  a  time-limited  function  whose 
Fourier  transform  best  approximates  a  bandlimited  function.  This 
approach  leads  to  the  well-known  Kaiser  window  expressed  by: 


Ig  (S/I  -  (K/N)2) 

r;T6r 


i(K)  = 


•N  <K  <N 


where  S  is  a  constant  (1  <  0  <  10) 

Iq  is  the  modified  Bessel  function  of  order  zero 
N  is  half  the  number  of  taps 

The  constant  0  can  be  determined  experimentally.  Figure  3-6  shows 
the  same  17  taps  frequency  response  with  Kaisier  window.  The 
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FIGURE  3-4.  LOWPASS  FILTER  RESPONSE  (17  Taps)  DIRECT  TRUNCATION 


FIGURE  3-5.  LOWPASS  FILTER  RESPONSE  (25  Taps)  DIRECTED  TRUNCATION 
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FIGURE  3-6.  17  -  Tap  LOWPASS  FILTER  WITH  WINDOW 


Pftssband  ripple  Is  essentially  vanished  and  the  stopband  rejec¬ 
tion  is  reduced  from  -20  dB  to  -38  dB,  Note  that  the  windowing 
operation  is  a  proper  scale  function  applied  to  the  impulse 
response  of  the  filter.  It  does  not  change  the  hardware  con¬ 
figuration  of  the  filter. 

3.4.2  A  Differentiator 

The  response  of  a  differentiator  can  be  written  by 
the  equation 

H(joj)  •  j(jj,  for  -ojj/2  <  (i)  <  0)^/2 

This  is  the  case  that  the  in-phase  response  is  zero  and  the  quad 
phase  is  an  odd  function.  As  a  result,  the  coefficient  (or  the 
impulse  response)  is  real  and  odd  functions;  that  is, 

c.  =  —  /  0)  sin  (ku)T)  dw 

s  0 

1  2_ 

=  /  oj  sin  (kwT)  du)  where  T*—  =Tr,u)  »2 
0  '^s  ^ 

“  -  ^  -1,1, ...N 

k  =  0 

where  k  *  0,  =  0. 

Figure  3-7  is  the  response  of  the  filter  with  17  taps  and  Figure 
3-8  is  the  response  with  Kaiser  window.  The  parameter  B  for  the 
Kaiser  window  is  set  to  3  to  reduce  the  ripple  due  to  direct 
truncation  of  the  Fourier  series  at  N  *  8. 
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IDEAL 


FIGURE  3-8.  (17  Tap)  DIFFERENTIATOR  RESPONSE  WITH  WINDOW 
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Hilbert  T ransformer 


The  frequency  response  of  a  Hilbert  transformer  can 
be  defined  as 


H(  juj) 


0  £u)  <  0)^/2 

-^  <  0.  <  o^s 


The  discrete  Fourier  transform  can  be  expressed  by  the  impulse 
response 


^  \  s i n'^ ( n7r/2 ) 


r\  f  0 

n  =  0 


The  comparison  between  a  17  tap  filter  and  the  ideal  transform  is 
shown  in  Figure  3-9. 

3.4,4  Equalizer  for  an  Arbitrary  Function 

The  transfer  function  of  an  arbitrary  frequency  response 
H(jw)  =  +  Q(ju)) 

over  the  band  of  interest  is  considered.  Since  the  samples  are 
not  symmetrical  in  any  sense,  the  resulting  tap  coefficients  will 
be  complex  and  asymmetrical.  The  matching  of  the  ideal  versus 
17  tap  equalizer  is  plotted  in  Figure  3-10  and  the  resulting  com¬ 
plex  taps  are  tabulated  in  Table  3-2.  Note  that  the  equalizer  taps 
are  complex: 

Ck  =  Xk  +  j  Yk*  k  =  1  ,2 . 17 
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RESPONSE  OF  A  (17  Tap)  HILBERT  TRANSFORMER  (EXAMPLE  3) 


and  the  gain  and  phase  adjustment  can  be  found  by  the  transforma¬ 
tion 

*  10  X  log  (X|^^  +  Y|^^) 

»  tan“^ 


TABLE  3-2.  RESULTING  EQUALIZER  TAPS 


Index 

Y  . 

1 

1 

0.438043 

-0.036516 

2 

0.931036 

0.043327 

3 

-0.092596 

-0.064665 

4 

0.086571 

-0.686610 

5 

-0.074379 

-0.578472 

6 

0.132042 

-0.643479 

7 

-0.213581 

-0.047679 

8 

3.133260 

0.643330 

9 

2.183238 

4.936495 

10 

3.120918 

0.712413 

11 

-0.187885 

-0.225917 

12 

0.090602 

0.764618 

13 

-0.032419 

0.508859 

14 

-0.005370 

0.732557 

15 

0.036480 

0.031661 

16 

-0.890202 

-0.018254 

17 

-0.468842 

0.016737 
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4.  CONTROL  PROGRAM 


A  Fast  Fourier  Transform  (FFT)  subroutine  is  employed 

to  calculate  the  set  tap  coefficients  for  the  equalizer.  The 

equalizer  taps  are  calculated  in  a  recursive  manner. 

The  input  to  the  AMTE  consists  of  a  set  of  N  complex 

correlations  as  a  function  of  frequency.  These  samples  can  be 

expressed  by  the  transfer  function 

+  j  Q(a).),  i  =  1  .2, - ,N 

t  h 

where  uj^  is  the  i  radian  frequency,  I(u)^)  is  the  in-phase 
frequency  response,  and  Q(a).)  is  the  quad-phase  frequency  response, 
and  N  is  the  total  number  of  sample  points. 

The  output  of  the  AMTE  is  a  set  of  equalizer  taps. 

For  the  system  at  hand  is  set  equal  to  17.  These  complex  tap 
coefficients  are  obtained  by  a  truncated  Fourier  series  coeffic¬ 
ients.  The  output  tap  coefficient  will  be  presented  by  its  mag¬ 
nitude  setting  (in  dB's)  and  phase  setting  in  degrees. 

4 . 1  FLOW  CHART 

A  simplified  flow  chart  of  the  main  routine  is  shown 
in  F i gure  4-1 . 

4.2  PROGRAM  IMPLEMENTATION 

The  program  is  implemented  on  the  RADC  HP  2100,  with 
the  operating  system  configured  on  2/6/81.  The  steps  to  be 
taken  for  execution  of  the  program  are; 

After  the  logs  on  the  HP  2100,  the  system  will  return 
the  prompt  sign  At  this  level,  it  is  necessary  to  link 

the  compiled  version  of  the  program  $ZA;:45.  To  link,  one  gives 
the  command: 
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FIGURE  4-1.  PROGRAM  FLOW  CKART 
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OF,  AIL  <CR> 


Then  after  computer  responds  with  then  type 

RU,  LOADR,  XX  <CR> 

where  XX  is  the  peripheral  number  of  the  user’s 
terminal  (this  may  be  optional)  and  <CR>  is 
the  carriage  return  key 

The  computer  responds  by  giving  a  prompt  of  “/LOADR:".  After  this 
prompt,  the  user  types: 

RE,  .  ZA  <CR> 

The  computer  then  lists  the  main  line  and  the  subroutines  and 
functions  called  for  by  the  user's  program,  in  this  case,  $ZA::45 
while  I’nking  them  together.  When  the  computer  is  through  linking 
these  together,  it  comes  back  with  another  ’’/LOADR:''.  The  user 
should  then  type 

/£  <CR> 

The  computer  than  links  and  lists  on  the  user’s  terminal  as  it  is 
linking,  library  routines  and  functions  necessary  for  executing 
the  user’s  program.  The  computer  then  comes  back  with: 

’’XX  PAGES  RELOCATED"  ’’XX  PAGES  REQ’D" 

"NO  PAGES  EMA"  "NO  PAGES  MSEC" 

’’/LOADR:XX  READY  AT  XX  (date,  time)’’ 

"LOADR-  $END  and  it  gives  a  prompt 
At  the  prompt,  the  user  should  run  the  program  that  has  been  loaded 
and  linked,  now  called  "AIL:.  The  command  is: 

"RU,  AIL"  <CR> 
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The  computer  then  outputs  anv  messages  and  data  from 
the  program,  as  it  is  running.  The  sequence  of  output  should 
be : 

INITIAL  SETTINGS  AT  STEP  0  SHOULD  BE... 

ENTER  THE  FILE  NAME 

At  this  point,  the  user  types  in  the  name  of  the  input 
file  (containing  a  distorted  sinusoid,  if  created  by  ATIO  as  des¬ 
cribed),  for  example  we  could  use  AD000,  one  input  file  from  ATIO; 
to  use  AD000  one  types  "AO000'‘  (return) 

after  the  message  "ENTER  THE  FILE  NAME". 

The  computer  then  prints  out  this  message; 

INTERMEDIATE  SETTINGS  AT  STEP  1  ARE: 

The  computer  then  prints  out  the  attenuator  and  the  line 
stretcher  values  for  the  equalizer,  then  continues. 

If  the  distortion  has  not  been  compensated  in  the  first 
iteration,  the  program  goes  back  to  subroutine  INPUT.  The  program 
repeats  this  loop  10  times,  or  until  the  distortion  is  less  than  a 
threshold  set  to  be  -45  dB,  whichever  comes  first.  When  the  pro¬ 
gram  reaches  the  point  of  exiting  the  loop,  it  prints  out  either: 

"AFTER  10  ITERATIONS  DISTORTION  CANNOT  BE  COMPENSATED" 


or : 


then ; 


"DISTORTION  COMPENSATED  AFTER  XX  ITERATIONS" 

"AIL:  STOP" 


The  computer  then  returns  to  the  system  level,  and  the  user  can 
perform  other  tasks,  or  sign  off. 


1 


A  sample  run  of  the  program  implemented  with  input 
data  is  contained  in  Appendix  1. 

4 . 3  COMPUTER  SUBROUTINES 

The  following  discussion  includes  the  entry  points, 
input  arguments  and  output  arguments  of  each  subroutine  in  the 
AMTE  program,  and  what  the  subroutines  or  functions  accomplish. 


1;  Subroutine  INTLZ  (DMYl,  DMY2,  N,  NTAP,  BETA,  IIO,  10,  KSW) 
initialises  constants  and  arrays. 

Input;  None 


2  .• 


n 

utput : 

DMYl  ,  DMY2 ,  M.  BETA,  10,  KSW 

DMYl 

A 

Initial  attenuator  (dB)  setting 

DMY2 

3 

Initial  phase  (degrees)  setting 

N 

Number  of  points  to  be  transformed 

NTAP 

Number  of  taps 

BETA 

Constant  used  in  KASER  subroutine 

10 

Output  device  number 

KSW 

Format  Flag 

KSW  =  0:  the  data  is  expressed  in  decimal  form 
KSW  -  1:  the  data  is  expressed  in  octal  form 

IIO  Input  device  number  (mag  tape  drive) 

Subroutine  DIAL  (NTAP,  DMYl,  DMY2,  X,  Y) 

This  subroutine  calculates  attenuator  and  line  stretcher  settings. 
Input:  NTAP,  DMYl  ,  0MY2 

DMYl,  DMY2  -  updated  attenuation  (dB)  and  phase  (degrees) 
Output:  X,  Y 

X  updated  line  stretcher  (division)  settings 

Y  updated  attenuation  (dB)  settings 
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I 


1 


1 


! 


1 


1 


I 

! 


3:  Subroutine  Head  (X.Y.NTAP,  ISW,  ICNT,  10) 

This  subroutine  displays  attenuator  and  line  stretcher  settings. 
Input:  X.  Y,  NTAP,  ISW,  ICNT,  10 

X,  Y  line  stretcher  and  attenuation  settings 

NTAP  number  of  taps 

ISW  ISW  =  0;  sidelobe  distortion  <  -45  dB 

ISW  =  1 ;  at  least  one  sidelobe  >  -45  dB 
ICNT  iteration  counter 

4.  Subroutine  INPUT  (N,  IBUF,  IIO,  10,  KSW) 

INPUT  defines  input  array  of  in-phase  and  quadrature-phase  data. 

Input:  N,  IIO,  10,  KSW,  file  name 

File  Name  file  on  which  the  data  is  stored 

Output:  IBUF 

IBUF  an  array  containing  both  I  &  0  data  in  octal  form 

5.  Subroutine  CHNG  (N,  IBUF) 

This  routine  swaps  the  upper  half  of  the  FFT  input  array  with 
the  lower  half. 

Input;  N,  IBUF 
Output:  IBUF 

6.  Subroutine  FORM  (N,  IBUF,  A) 

FORM  will  compute  complex  array  A  =  I  +  jQ  and  also  gain  = 

5  log.|Q  (I^+Q^)  and  phase  =  tan”^  (Q/I) 

Input:  N,  IBUF 

Output:  A 

7.  Subroutine  INVT  (N,A) 

This  subroutine  computes  the  sequence  to  be  equalized: 
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Input: 

Output 


A .  =  i—  where  i  =  1  ,2 , 3  .....  N 


N  ,A 


10, 


1 1 


Subroutine  FFT2C  (A.M.IMK)  computes  the  Fourier  transform  of 

the  input  array  A 

Input:  M , A 

Output:  A 

Subroutine  TAPCO  (N.  A.  NTAP,  COEF,  10) 

TAPCO  truncates  array  A(N)’  into  finite  terms.  The  number  of 
terms  is  defined  as  NTAP  (number  of  taps)  to  specify  an  equal¬ 
izer.  The  variable  10  is  a  device  number  for  the  output  which 
is  normally  computer-dependent.  The  resulting  equalizer  taps 
are  stored  in  COEF  (NTAP). 

Subroutine  KASER  (NTAP,  BETA,  COEF,  10) 

This  subroutine  applies  the  KAISER  window  function  to  the  trunc¬ 
ated  Fourier  transform  output  COEF  (NTAP).  A  functional  sub¬ 
routine  BESSL  (X)  is  required  to  compute  the  window  parameters. 
BESSL  is  the  modified  Bessel  function  of  the  first  kind  with 
zero  order.  The  argument  for  the  BESSL  is  X. 

Subroutine  OUTPT  (DMYl,  DMY2,  COEF,  NTAP,  ISW) 

Subroutine  OUTPT  computes  updated  attenuation  (dB)  and  phase 
shift  (degrees) . 

X.  =  20  log^g  do'');  A  =  +  1 


Y.  =  20  log^Q  do”);  B  = 


20 

ATT. 

20 


+  1 


updated  attenuation.  =  20  log,^  (10^) 
X.+Y.  ^  ''' 

where  C  =  -g-Q  ^  +  1 

and  i  =  1,2,....N 
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DMYl  -  previous  attenuation  (dB) 

ATT  -  present  tap  indication  (dB) 

phase  shift  = 

i=l  ’ 

Input:  DMYl.  DMY2,  C0£F.  NTAP 

DMYl  -  updated  attenuation  setting  (dB) 

DMY2  -  updated  phase  shift  setting  (degrees) 

ISW  -  ISW  =  0;  si  delobe  distortion  <  -45  dB 

ISW  =  1;  at  least  one  sidelobe  >  -45  dB 

12:  Subroutine  DIAL  (NTAP,  DMYl,  DMY2,  ISW,  ICNT,  10) 

DIAL  was  described  previously. 

13:  Subroutine  HEAD  (X,  Y,  NTAP.  ISW,  ICNT,  10) 

HEAD  was  described  previously. 

14:  Subroutine  END  (ISW,  ICNT,  10)' 

This  subroutine  will  print  a  message  depending  on  what  the 
inputs  are: 

Input:  ISW,  ICNT,  10 

Output:  If  ISW  =  0,  message  is  "Distortion  compensated  after 

so  many  iterations" 

If  ICNT  =  10,  message  is  "After  10  iterations  distor¬ 
tion  cannot  be  compensated" 

15:  Function  ZL0G2:  This  function  calculates  log  to  the  base  10 

of  the  input  quantity  X.  It  takes  the  natural  log  (ALOG)  of  X 
and  multiplies  by  Y,  where  Y  =  l/log^lO. 

16:  Function  ZNl :  This  calculates  arctangent  of  input  quantity 

(in  radians). 
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4.4 


COMPUTER  PROGRAM  LI’^TING 


The  following  is  a  listing  of  the  computer  program; 


0  0  0 1 
0  Cl  Oc 
000? 
00  04 
0005 
0  0  Cit- 

00  or 

ij  0  ij? 
000? 
ii  Ci  1  0 
0  0 1 
0  C- 1  ? 
0  C.  1  J 
0  0 1  4 
U  Ij  1  5 
Cl  0 1  6 

0  0 1  r 
0018 
0  0 1 9 
0  08  0 
0081 
0  08S 
008? 
0084 
0  085 
0  08b 
008? 
0028 
0029 
0  03  0 
0  031 
0038 
0033 
0034 
0035 
00  3  b 
0  037 
ij038 
0039 
0040 
0041 
0048 
004  ? 
0  044 
0045 


FTN4 

C  PEVICEH  FILE  1-88.'81 

PPDGFfiM  HIL 

C  aMFLE:-:  h  ■  858  >  .  E  <  85b:>  •  CDEF  'i.85b> 

Ii  I  mens  I  on  I  EUF  •:  85b.>  •  DMV 1  C 1  7  >  .  Df1Y8  C 1  7>  .  X  (.1 7>  »  Y  <.  1  7> 

DIMEMSinri  IWf  19> 

ICNT=0 

M=8 

C  IMITIhLIZE  COMSTRNT'- 

CftLL  inTLZ  (IiMYl.IiMYS.M.hTHF'.EETH!.  IIQ.  ID.K8W,' 

C  CftLCULRTE  IhTITIRL  hTTENUhTDF  RMH  LINE  STFETCHEP  SETTINb 

C  RLL  D I RL 'NT RP «  DNY 1 • DMY8  *  X »  Y) 

C  ni'ZFLRY  RTTENURTDP  RND  LINE  STRETCHER  SETTINGS 

C  RLL  HERD  ■: : ;  *  Y »  NTRP »  I SW ■■  I CNT »  I  □> 

DD  10  I =  1 »  1  0 
!i:nt=  I 

I‘^<ICNT.EQ.8;'  GO  TD  80 

C  ENTER  SRNPLE  TRRNSFER  FUNCTION  RS  I  RND  Q  DRTR 

CRLL  INPUT  <N«  I £iUF ?  IIO»  10»KSW  ‘ 

C  SHIFT  DRTR  N  8  SRtlPLES 

CRLL  CHNGcN. IEUF> 

C  COMPUTE  TRRNSFER  FUNCTION  R  CN.'.' =I  CNi +JQ  (N;' 

CRLL  FORM  cN.  IEUF-R;- 

C  COMPUTE  EQURLIZER  TRRNSFER'  FUNCTION  l-'R^.N-' 

CRLL  INVT(Ni.R> 

C  PERFORM  FFT  ON  R 

CRLL  FFT8C < R ’ M » I Wk > 

C  TRUNCRTE  THE  IMPULSE  RESPONSE  TO  17  TERMS 

CRLL  TRPCO  'N.R.NTRP^COEF.  ID;' 
r  PERFORM  I'JINDOW  FUNCTION  ON  TRUNCATED  SERIES 

CRLL  KRSER  ■ NTRp. EETR- COEF> 

C  CRLCULRTE  RTTENURTION  <Dt:;'  RND  PHRSE  SHIFT  (DEGREES;' 

CRLL  OUTPT  '  DMYl  ■.DMY8»CDEFjNTRP- ISW.' 

C  CRLCULRTE  RTTENURTION  <TURNS>  RND  LINE  STRETCHER  SETTINGS 

CALL  DIAL  <  NTRF>  DMY 1  •  DMY8'*  Y:' 

C  DISPLAY  RTTENURTION  RND  LINE  STRETCHER  SETTINGS 

CALL  HERD  (K* NTRP»  ISW?  ICNTj  IO> 

IF'-ISW.EQ.  0;'GO  TO  80 
10  CONTINUE 

C  DISPLAY  AFTER  HDM  MANY  ITERATIONS  DISTORTION  I'lRS 

C  COMPENSATED  .DR  THAT  IT  MRS  NOT  COMPENSATED  AFTER  A 

C  GIVEN  NUMBER  OF  ITERATIONS. 

80  CRLL  END  '  I  SM.  ICNT*  lO;- 
9998  I'lR  I  TE  '37. 9998 0  > 

99980  FORMAT'"  AIL  s  STDP"> 
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0047 
0048 
0049 
0050 
0051 
0058 
0058 
0054 
0055 
0058 
0057 
005S 
0059 
0080 
0081 
0088 
0083 
0084 
0085 
0088 
0087 
0  088 
0089 
0070 
0071 
0078 
0073 
0074 
0075 
0078 
0077 
0078 
0079 
0  08  0 
008 1 
0088 
0083 
0084 
0085 


C  CfiLCUftTE  LOG  TD  THE  EftSE  10 

EUNCTioH  ZLOGi  <>;;■ 

V= 0.43489448 19 
IFCX.LE. 0.>K=1.E-10 
ZLDGl=Y*i=(LOG  OO 
FETUPN 
END 
C 

C  CflLCULPiTE  HRCTFINGENT 

C 

FUHCTIDN  ZN8<:Y'.X> 

I F  <  X . EQ . 0 . >  X= 1 . E- 1 0 
FI =3. 141598854 
Z=YVX 

ZN8=FlTFiN  '.  Z'' 

I F  <:X .  LT .  O'-  ZNE'^ZHa+P I 
C  WFITE  <37..978>  Z..ZN8 

978  FDFM8T  a  X, E18. 5» 1X» E18. 5) 

FETUFN 

EMD 

C 

c 

C  CFLCULftTE  EESSEL  FUMCTIDM  DF  THE  FIRST  KIND  WITH  0  ORDER 

FUNCTION  EESSL<X> 

Y=X^'8 

DELTft»lE-8 

E=l. 

DE=1. 

DO  1  1=1.85 
DE=DE^Y-'FLOHT<:i> 

SDE=DE*IiE 

E=E+S:r'E 

IF<E*DELTH.GT.SDE>GOTO  10 
1  CONTINUE 
1 0  EESSL*E 
RETURN 
END 
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SWAP  THE  UPPER  HALF  OF  THE  ARRAY  WITH  THE  LOWER  HALF 
SUBROUTINE  CHNG CN? IBUF> 

INTEGER  IBUF  <1  >  » IBUFl  <:300> 

K.l  =  <N+l>/a 
K2=K1-1 
K:=N."c' 

DO  10  i=k:i,n 
J=I-K 

IBUFl <J>=IBUF(I) 

CONTINUE 
DO  £0  I»1.K 
J*I-t>Kl 

IBUFl  <J)=IBUF<:i> 

CONTINUE 
DO  30  1  =  1  »N 
IBUFa)=IBUFl  <I> 

CONTINUE 

FORMAT  <lXf l£<I6t 1X>) 

ZX=1 

CALL  TYPE<ZX) 

RETURN 

END 


! 


0110 

C 

01 1 1 

c 

01 13 

c 

0113 

C: 

0114 

c 

01 15 

0116 

0117 

0 1 1 8 

01 19 

0130 

0 1  £  1 

0133 

0133 

0134 

0135 

0136 

0137 

0 1 38 

0139 

0130 

0131 

0133 

0133 

0134 

0 1 35 

0136 

0137 

0133 

0139 

0140 

0141 

0143 

1  0 

0143 

c 

0144 

C 

0145 

979 

0146 

0147 

0148 

0149 

PEPFDPh  POLYNDMIhL  FUNCTION  TD  CONVEPT  FPOM  hTTENUmTIDN  *DE:' 
TO  hTTENUhTION  '  DIhL  TUPNS>  FtND  FPDM  PHFiSE  SHIFT ':DE6PEES  > 
TD  LINE  STPETCHEP  SETTINGS  dilVISIDNS;' 

SUEPOUT I NE  n  I  FiL  CNTFiP ,  ErMY  1 .  DMY£  ■■  X •  Y> 

D I  MENS:  I DN  H  i.i '  1 7 >  R 1  <.  1  T  j  »  fic.’  C 1 7>  •  RS  ( 1 7^‘  • 

1  DMY 1  <  1  >  I.  riMY£  <  1  >  .  X  a  >  »  Y  <  1  > 

DRT h  RO  '  “ .  '?•  1  OS' •  — .  7£b  1  •  — .  ■r'c'5 1  >  — .  GSS'  1  >  ~ .  8t  07  !>  —  1 .  £S"?7 « 

1~1.  Icl6»~l.  Ic'lG  .  y644 »  — .  £:4 1  ~  1 .  OS'S  1 »  — .  ■?  1 1-9  > 

S'-S.  1  4  1  3 •  “ .  84c'c' ’  -  1  •  •j:34c'f  7S;43:x 
DRTR  R 1  .  1  1 66 »  - .  07'3c'!ri »  1 1 1 6»  ~ .  06SSS »  ~ .  1  0S6 »  “ .  1  S'SSS » 

1-.  1  07Si«  1  075?  0.  OSS 05?-.  07806?  -  .  1176?-.  0'S46?-.  IS'S'S? 

S'- .  05 1  78  ?  - .  1 1  '?S'  ?  - .  0777 1 

DRTR  RS-'-.  0009174? .  OOIS'95?  0003383?  .  001839?  .  0003631  ? 

1 . 0  0  0  071 6  5?-.  0  0  03634 ?  - .  0 0 03634 ?  0 .  ? .  0  0 1>3  096  ?  - .  0  0  03863  ? 

£-.  00003543?  .  000673? .  0006743?  .  003387?  -.  001 177?  .  OOIS'85'- 
DRTR  R  3  .  0  0  0  0 3 8 3797?.  0  0  0  0  0377?-.  0  0  0  0 3436?.  0  0  0  0 1 1  0  7 5  ? 

1-.  0000313?  -.  000033518?  -,,00003773?  -.  00003775?  0,  ?  -.  0O00S50'56? 
£-.  00003957?-.  0000373?-.  0000166'96? -.  00001571?  .  00001376? 

3-.  00004344?-.  00000'3316^' 

DO  10  I=1?NTRP 
J=I 

IF<I.LE.8>  J=9-I 

Yl=R0<:i>+Rl  i:l>*DMYl  <:.J> +R3  <:i>  ♦DMYl  (J>«'DMY1  (J^ 

Y£=R3  CD  ♦DMYl  CJ;-  ♦DMYI  <  J>  ♦DMYl  U,> 

T  <.I>  =  V1+Y3 

X  <  I>  =DriY3  ( J>  ♦5.634-'  i360.  ♦!  .  09 1> 

IF  CY  <I>  .  LT’.  0.  0>  Y  <  I>  =  0.  0 
I F  <  Y  <  I  >  .  GT .  1  0 .  0>  Y  <  I  >  =  1  0 .  0 
IF  <:X  i:.  I>  .  LT.  0.  0>  X  <I>  =0.  0 
IF  cxa>.GT.  11.  o:-  xa>=i  1.0 
CONTINUE 

WP I TE  <  37  ?  979>  <;X  v  I  >  ?  I  *  1  ?  NTRP) 

UP  I TE  <  37?'579^  <  V  <  I  >  ?  1  =  1  ?  NTRP) 

FDPMRT  <1X?4<E1£.5? 1X?E13.5? 1X>> 

ZX=£ 

CRLL  TYPE  (;ZX> 

PETLIPN 

END 
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0 1 5  0  C 
i;i  151  c 

015£  C  PRINT  DLIT  WHETHER  DISTORTION  IS  COMPENSRTED  ftFTER  THE 

015:-:  C  LfiST  ITERHTIDN 

0154  SUDRDUTINE  END  <  I  :  W.  IC:NT»  IO> 

0 155  I F  <;  I SW .  EO .  0>  GO  TO  10 

0 1 56  WP I TE  •:  3 ' >  3  0;'  I NC T 

0157  £0  FORMHT<.',3>:.5HFtFTER.  I3*31H  ITERATIONS  DISTORTION  CAN  NOT. 

015S  114HEE  COMPENSATED*-  '' 

0159  GO  TO  30 

0160  10  CONTINUE 

0161  ICNT=ICNT-1 

0 1 62  WR I TE  <  37 • 4  0> I CNT 

0163  40  F0RMAT<.*3X*£SHDIST0RTIDN  COMPENSATED  AFTER* IE> 1 IHITERAT ION 

0164  30  CONTINUE 

0165  ZX*3 

0166  CALL  TYPECZX;' 

0167  RETURN 

0168  END 

0169  C 
0170  C 
0171  C 

017£  SUBROUTINE  FORM  <N!.  IBUF.  A> 

0173  COMPLEX  A<1> 

0174  DIMENSION  II:UF<1> 

0175  C  COMPUTE  COMPLEX  ARRAY  AsI-*JQ 

0176  DO  10  1=1 >N 

0177  IA1=IEUF  <  i;- --£56 

0178  IE1=IEUF<I>-£56»IA1 

0179  IFflAl.GE. l£S> IAl=IAl-£56 

0180  IF<IE1.GE. IEl=IBl-£56 

0181  XI=FLOAT  <  IA1> -- 128. 

0182  XQ=FLOAT  <  lEi:- -"1£8. 

0183  A  <  I  >  =CMPLX  <  X I  >  XG!) 

0184  10  CONTINUE 

0185  998  FORMAT  •  IX*  4  <E1£.  5*  1X»E1£.  5-.  IX)  > 

0186  ZX=4 

0187  CALL  TYPE<ZX) 

0188  RETURN 

0189  END 

t.l«'0  C 
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T 


0191 

0 1 9£ 

c 

3LIBRDUTIME  HEFiD  V.  NTfiP»  I SWi.  IC NT. 

0 1 93 

D I  MENS  1  ON  X  1 1  .  Y  a  .  ZR  a  7> 

0194 

DhTh  ZR'4H1  h=.4H3  H=>4H:3  R=<4H4  H=.4H5  M=j4H6  R=!i4H7  R=> 

019^ 

14H8  R=.4H  M  =.4H1  B=.4H3  E=.4H:3  E=!.4H4  E.=»4H5  B=.4H8  £=• 

0 1 9t. 

14H7  E='.4HS  E=' 

0197 

I.IPITE  t  SZ.  145  :' 

0 1 93 

145 

FDRNRT ' IHl > 

0 1 99 

Ii:NTl  =  Ii:NT-t-l 

1 !  3  i;i  0 

i.Q  TO  '  1  0 . 3 0 .  c! 0 . 3 i.i » 3 0 .  c 0 » 3 0 » 3 U j  3 1.1 . 3 0 »  Z"  1.0’  ICNTl 

Oc'Ol 

1 1'l 

CONTINUE 

0303 

I’lRITE  '  37.  150:'  ICNT 

Cl  3  03 

1* 

DISPLRY  RTTENURTION  RND  LINE  STRETCHER  SETTINGS 

0304 

1 5  0 

FORMPTc  '-.  18X.35HINITIRL  SETTINGS  RT  STEP  .Il.llH  SHOULD  EE: 

03  05 

GO  TO  305 

03  03 

30 

COr<TINlJE 

0307 

I  Sl.l  1  =  I  St’i+ 1 

0308 

gD  TO  '■  i.i .  c'  i.i  i.i  ■'  I S  l-'.i  1 

03  Cl  9 

3  0  0 

CONTINUE 

03 1  0 

WRITE  <  37.  155> ICNT 

03 1 1 

155 

FORMRT’O'.  18X.39HINTERMEriIRTE  SETTINGS  RT  STEP.I3»5H  RRE:  > 

0313 

GO  TO  305 

03 1  3 

3  0 

CONTINUE 

03 1  4 

WRITE  <37. 18 CO'  ICNT 

03 1 5 

1 8  0 

FDRMRT<.'.31X!>33HFINRL  SETTINGS  RT  STEP  »I3>5H  RRE:  > 

03 1 8 

305 

CONTINUE 

03 1  7 

i.iRITE<37.985:’ 

03 1 8 

985 

FDRMRT  1  •  1  OX.  1  7HRTTENURTOR  <TURNS>  »31X.  19HLINE  STR£TCHER  i^DI  V 

03 1  9 

DO  310  I=1.NTRP 

0330 

lilR  I TE  <  37 .  17 CO'  ZR  <  I  /  » Y  <  I  >  «  ZR  <  I  ^  *  X  <. I ) 

033 1 

170 

FDRMRT  <  1  OX .  R4.  F8 . 3.  3i4Xl!i  R4!i  F4.  1  > 

0333 

310 

CONTINUE 

0338 

ZX=5 

0334 

CRLL  TYPE  <ZX) 

0335 

RETURN 

0338 

0337 

C 

END 
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Ci££8  C 

OSSS*  SUERDUTIMF  INPUT  CN,  lEUP..  1 10»  10.  KSW> 

Oc3  ti  n  I  MENS  I  ON  I  nCB  •:  £5t.>  .  NRM  -  SO  »  I BUF  ( £5fc.> 

t3£3 1  D I  MENS  I  ON  I  BUF  0  <  1  £8>  .  I  BUF  1  <  1  £S> 

0£3£  WRITE  <3?.  1  O:' 

0£33  10  FORMPTc/-  .  '  ENTER  THE  FILE  NHME  '  > 

0£34  REHD(37.  15>Nl=tM 

0£35  IS  F0RMFiT<3H£) 

0£3F  IL=M 

0£37  CftLL  □PEN<  IDCB.  IERR.NFtM> 

0£38  IFdERR.LT.  0>i50  TO  900 

0£39  C  ENTER  BHTft  IN  OCTHL  FORMFiT 

0£40  CRLL  RERDFCIDCB.  lERR.  IBUFO;> 

0£4 1  CRLL  RERIiF  <  I BCB .  I  ERR .  I  BUF  1  > 

0£4£  IFdERR.LT.  0>  GO  TO  910, 

0£43  CRLL  CLOSE d BCB. lERR) 

0£44  C  fill  RRRRY  WITH  IN-PHRSE  RND  GiURDRRTURE-PHRSE  COMPONENT 

0£45  DO  60  J=l.l£8 

0£46  I BUF  <  J  ' = I EUFO  C J> 

0£47  60  CONTINUE 

0£48  DO  70  K=l,l£8 

0£49  lEUF <K>1£8>=IBUF1 CK) 

0£50  70  C:DNTir<UE 

0£51  998  FDRMRT  dX.l£d6.1X)> 

0£5£  GO  TO  999 

0£53  900  WRITE <37. 30> lERR 

0£54  30  FDRMRT  C-'  .  "  FMP  ERROR  "14) 

0£55  GO  TO  999 

0£56  910  WRITE  <37.  £0:> 

0£57  £0  FDRMRTf'.'  ERROR  RERDING  THE  FILE  ') 

0£5S  999  ZX=6 

0£S9  CRLL  TYPE  CZX) 

0£60  RETURN 

0£61  END 

0£6£  C 
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0£63  C 
i:i£64  C 
Ocb-' 

Oce-b 
02  6  7 
0£b8 
0269 
027  0 
0271 
0272 
0273 
0274 
0275 
0276 
0277 

0278  1 0 

0279 
0280 
0281 
0282  C 
0283  C 
0284 
0285 
0286 
0287 
0288 
0289  C 


INITIRLIZE  COMSTRNTS  HND  RPRRY 
S  UEROUT I  ME  I MTLZ  <  DMY 1 »  DMY2 .  h.  HTRP.  BETR .  1 1 0»  K:SW> 

D I  riENi;  I  OH  DMY 1  <  1  ;:•  »  DMY2  C 1  >  » R  Cl  ?>  B  c  1 7:> 

DRTR  R-  -59. 964 • -59. 964 « -51 . 072? -56 . 743 ? -77. 661 > -6  0. 55 1 • 
1  —  72 .  7 1 4  ?  —  57 .  5o  ?  0 .  ?  —  52‘ •  449 ?  — 6c‘ .  1 4-^'  ?  —  55 .  3 1'  3  ■*  —  5  .  r  4 1  ? 

2-54 .419?  -69 . 5  04  ?  —57 .  047  ?  —58 . 965-' 

DRTR  Ei'''659 . 23  ?  756 . 49?  2: 1  0 . 52i  ?  982i .  44  ?  162.  1  ?  399 . 86  • 

1670.  ij3'*2'59.  16?443.  U9?  129,  6C'»  399.  S6’*  561 . 96'?  778.  1? 

254.  04?  151 . 3?  345.  S2?670.  03/'' 

N=256 

NTRP=17 

DO  10  1=1?NTRP 
DM VI  Cl)=Ra> 

Df1Y2  d>  =BCI> 

CONTINUE 
EETR=1. 0 
1 10=37 
10=37 

FOR  DECIMRL  I«,G!?KSUI=0 
FDR  DCTRL  IS:Q?KSW=1 
KSW=1 
ZX=7 

CRLL  TYPECZX> 

RETURN 

END 
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1 


0c'90 

C 

0291 

C 

CftLCULftTE  THE  E'ECIPPOCftL  DF  THE  TFHHSFER  FUHCTIDN 

0292 

SUBROUTINE  INVT<N.ft> 

0292 

COMPLEX  l=(a> 

0£94 

DO  10  I=1.N 

0295 

ft  <  I  >  =  1 .  X  R  U  :•  ♦FLOftT  c N> 

ijc'9b 

1  0 

CONTINUE 

i:i£97 

c 

IJR I  TE  C37 «  998)  <  ft  <I)»I  =  1»  N> 

0£98 

FOPMftT  ^  1X!.4<E1£.5!.  1X.E1£.5.  1X>> 

ij£99 

ZX=8 

0300 

CftLL  TVPE<ZX) 

0301 

PETUPN 

03  0£ 

END 

0303 

c 

0304 

c 

0305 

c 

ftPPLY  KftlSER  MINDOW  FUNCTION  TO  SMOOTH  OUT  THE  RIPPLE 

03  06 

c 

TO  TRUNCftTION  OF  THE  FOURIER  SERIES 

0307 

SU£;RDUTINE  KftSER  <.N»B»C> 

0308 

COMPLEX  CCD 

0309 

IF  <  B.  EQ.  0.  0)130  TO  9999 

0310 

NV£=<N-*-1)'^£ 

031 1 

DO  10  I*DN 

031£ 

Z£=FLDftT  CI-NV£)  ^  FLOftT  CN.'£) 

0313 

Z3=B*SG'RT  Cl .  -Z2*Z2'j 

0314 

CCD  =C  C I )  ♦BESSL  CZ3)  EESSL  <.B> 

0315 

1  0 

CONTINUE 

0316 

c 

WRITE  C37.979)  CC  CD  I  =  DN) 

03 1 7 

979 

FORMftT  <  1X.4C1X!.E1£.55  1X*E1£.5)) 

0318 

9999 

ZX=9 

0319 

CftLL  TYPECZX) 

0320 

RETURN 

0321 

END 

0322 

C 
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0?c3 

i:i3c4 

0335 

0336 

0337 

0333 

0339 

033  0 

0331 

0333 

0333 

0334 

0335 

03  36 

03:37 

0338 

0  339 

034  0 

0341 

034  3; 

034  } 

0344 

0345 

0346 

0347 

0348 

0349 

0350 

035 1 

0352 

0353 

0354 

0355 

0356 

0357 

035S 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0368 

0369 

0370 

0371 

0372 

0373 


C 

:  UBP'DUT  I  hE  OUTPT  <  DMY 1 .  OMYa.  C  OEF  .  NT8P .  I  SU;' 

COMPLEX  COEF-  l,' 

D  I  MEM  :  I  ON  DMY  1  <  !  >  .  DMY2  <  1  >  •  P(TT  a  7>  ,  PHS  i  1 7> 

DO  10  I=1.MThP 
XP=PEhL  >  CDEF • I >  > 

XI=8IMHi3  '  COEF  <  1  :•  - 

C  compute  hTTEMUhTIDM  of  DISTORTION 

ftR=XR^XR+XI^XI 
RTT  i  I>  =1 0. ♦ZLC61  <HH> 

C  COMPUTE  THE  PHhSE  SHIFT  OF  ThP  CX)EFF  IC  I  ENTS 

PHS  (' I ■>  =  .:n2  > i .  xr>  ♦45.  -- hthn  a .  > 

3Z=ZN8<>:i  .XR- 

L  I'lP'ITE  *j7«97c'':'  PHS<I>«ZZ 

978  FORMhT  ■  IX- E12.  5.  1X-E12.5:- 

10  continue 

DO  20  I=1.NTHP 
FiTT  I  '  =hTT  '  I  ■>  -hTT  C9> 

PHS  '  I  >  =PHS  1  •  -PHS  <9> 
cO  CONTINUE 
I  S  l,ii  =  0 

NTftPl=NTRP-l 
DO  30  I=1»NTRP1 
J=I 

IFa.6E.  9>  J=I  +  1 
IF<RTT(.J.' .LE.-45.>GD  TO  15 
I  Si.l=  1 
GO  TO  30 
15  CONTINUE 
ISU=ISU 
30  CONTINUE 

I F  I  I  S  l•J .  E Q .  0  ■'  GO  TO  40 
DO  50  I=1-NThP 

X=20.*ZLOGl  <  dO.^^'DMYl  <I>  20.;'>-*-i;' 

Y=2  0  ♦ZLOG 1  <  1  1  0 .  ♦♦  <  FiTT  ■,  1 0  -'2 0 .  >  >  + 1  > 

DMYi  ( 1  ■■  '20.  ♦ZLOGl  ■:  '  1  0.  ♦♦  <  (X!+Y-'  ••  20.  -l^' 

DMY2  <:  I  >  =DMY2  C  I  >  +PHS  <  I  > 

50  CONTINUE 
40  CONTINUE 

URITE  <:37:.996> 

996  FORMHT  ( 1 X. 1 4HHTTENUHTDR cDEO  > 37X, 1 9HL I NE  STRETCHER <DEG> . -u 
*  888  ZX= 1 0 

DO  SSI  I  =  1'.NTHP 

UR  I TE  <37 . 977>  DMY  1  <  I  >  .  CMYc'  <  I  > 

977  FORMHT  < 1 X . E 1 2. 5. 37X. El  2. 5> 

SSI  CONTINUE 

979  FORMHT  < 1 X . 4 v 1 X» E 1 2. 5* IX- E 12. 5> > 

CHLL  TYPE<ZX> 

RETURN 

END 

C 
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0374  C 

0375  SUEPDUTIME  TFiPCD  (HfR.NTFtPjCDEF*  !□> 

0376  COMPLEX  M  U  > - CDEF  < 1 > 

0377  NV3=  (;NTRP+1>  -•  3 

037S  NV21=M+NVc+l 

0379  NV1=MV3+1 

0330  C  TPUNCRTE  THE  RRRfiY  ft  CM)  INTO  R  GIVEN  NUMBER <NTRP>  OF 

0331  C  TERMS 

0333  DO  30  I  =  1:»NV3 

0333  J=NV3-I+1 

0334  COEFCD^RCJ) 

0335  £0  CONTINUE 

0336  C  OBTRIN  THE  UPPER  HRLF  OF  THE  COEFFICIENTS  CENTERED 

0337  C  RROUND  <NTRP+l>/'3 

0388  DO  30  I=NVl!.NTRP 

0389  J=NV31-I 

0390  CDEF<:i>=RCJ> 

0391  30  CONTINUE 

0393  C  WRITE  <37»979>  CCOEF  <!>  »  1  =  1  >  tHTRP> 

0393  979  FORMRT  <1X> 4 ( 1 Xj E13. 5» 1X» E13. 5> > 

0394  2X=11 

0395  CRLL  TYPE  C2X> 

0396  RETURN 

0397  END 
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0398 

0399 

0400 

0401 

04  03 

0403 

04  04 

04  05 

04  06 

0407 

0408 

0409 

041  0 

0411 

04 1  £ 

0413 

0414 

04 1 5 

04 1 6 

04 1 7 

04 1 8 

0419 

0430 

04£  1 

04  ££ 

04  £3 

04  £4 

04£5 

04£6 

04  £7 

0488 

04£9 

0430 

0431 

043£ 

0433 

0434 

0435 

0436 

0437 

0438 

0439 

0440 

0441 

04 4 £ 

0443 

0444 


C 

C 

SUBROUTINE  FFT£C 
INTEGER  M .  I WK  <  1  ;■ 

COMPLEX  Pa;* 

C 

C 

INTEGER  I  -  I SP .  J .  J J JSP .  K  » K 0 •  k  1  >  K £ .  K3 >  K.B ,  KN » MK  »  MM «  MP .  N » 
1  N4»N8'.N£.LM.NN.  Jk 

RE PL  RPD • C 1 • C£ ?  C o • S 1 ?  S£ »  £3 • Ck  >  S k  »  SQ  j  PO  j  P 1 • P£ >  P3 > 

1  B  0  <  B 1 . B£ - B 3 - TWOP 1 . TEMP . 

C  £  ZERO »  ONE  •  Z 0  <■  £’>  ?  Z 1  c£ •'  *  Z3  <  £^ 

£  ZERO ' ONE • Z  0  <  £>  »  Z 1 <  £ >  «  Z3  <£> • Z£  » £ > 

COMPLEX  ZPO?  ZPl >  ZP£«  ZP3-  Pk£ 

EQUIVPLENCE  (;ZP0«Z0  <  1>  >  *  (ZPl  •  Z1  <  1>  >  >  <:ZP£«Z£  (.r.>  >  . 

1  <ZP3 •  Z3  <  1  >  >  »  i. PO'  Z 0  <  1  >  >  »  <  E O?  ZO  '  £>  )  >  < P 1  >  Z1  1 1  v  >  > 

£  <  B 1  *  Z 1  ( £  j  /'  ?  (Pc'?  Z£  ( 1  >  >  •  (.B£ »  Z£  ( £)  ^  »  (.P3 ?  Z3  (.  1  -O  ? 

3  (B3  *33  (!£;•> 

DhTP  SQ.  .  7071  06781-' 

BPTP  S  k  .  38£6S343!£/'' 

DmTP  Ck .  '?£3;87953!3'^' 

DP  TP  TMDP I  6 .  £83 1 853 1 
DPT P  ZERO/  0  •  Cl/  >  ONE/ 1 .  Ci/ 

MP  *  M+1 

H  = 

lUlKCl;  «  1 
MM=<M/£>^£ 
kN  =  N+1 
DO  5  I=£?MP 

I  w  k  <;  I  >  =  iw  K  <  I  - 1  >  + 1  ui  k  c  I  - 1  > 

5  CONTINUE 

RPD  =  TUOPI/N 
Mk  =  M  -  4 
kB  =  1 

IF  CMM  .EG.  M>  GO  TO  15 
k£  =  kN 

kO  =  IWk  CMM+r.'  +  kB 
10  k£  =  k£  -  1 
kO  =  kO  ~  1 
Pk£  =  P<k£> 

P<.k£>  =  ft(K0>  -  Pk2 
PCKO)  =  P(kO)  +  Pk£ 

IF  CkO  .GT.  kB>  GO  TO  10 
15  Cl=ONE 

SI  =  ZERO 
JJ  =  0 
k  =  MM  -  1 
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0445 

J  =  4 

0446 

IF  CK  .GE.  1>  GO  TO  30 

0447 

GD  TD  70 

0448 

£0 

IF  <:iWk:<J)  .GT.  JJ>  GD  TO  £5 

0449 

JJ  =  JJ  -  iwk;(j> 

0450 

J  =  J  -  1 

0451 

IF  .GT.  JJ>  GDTD 

£5 

045£ 

J  j=jj-iwk:  (  j> 

0453 

J=J-1 

0454 

K  =  K .  +  £ 

0455 

GQ  TD  £0 

0456 

£5 

JJ  =  IUIK<.J>  +  JJ 

0457 

J  =  4 

0458 

30 

ISP  =  iwk<k:> 

0459 

IF  (JJ  .EQ.  0>  GO  TD  40 

0460 

C£  =  JJ  ♦  ISP  ♦  RftD 

0461 

Cl  =  CDS<C£> 

0463 

SI  =  SIM(;C;£> 

046? 

35 

C£  =  Cl  ♦  Cl  ~S1  ♦  SI 

0464 

S£  *  Cl  ♦CSl  +S1> 

0465 

63  =  C2  ♦  Cl  -S2  ♦SI 

0466 

S3*  C2  ♦SI  ♦  S2  ♦Cl 

0467 

4  0 

JSP  =  ISP  ♦  KB 

0468 

C 

WRITE  <;37»997>  Sl»S2fS3 

0469 

C 

WRITE  <;37!>997>  C1»C£»C3 

047  0 

DD  50  I  =  1  »  ISP 

0471 

KO  =  JSP  -  1 

0472 

K1  =  KO  ♦  ISP 

0473 

K2  =  K1  +  ISP 

0474 

k;3  =  K2  +  ISP 

0475 

ZftO  *  RCKO) 

0476 

ZFll  =  fi<Kl) 

0477 

Zfl£  *  fiCKE) 

0478 

ZFl3  =  fi<:k3> 

0479 

IF  (SI  .EQ.  ZERD)  GD  TD 

45 

0480 

TEMP  «  ftl 
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0481 

81  =  fll  ♦  Cl  -B1  ♦SI 

0482 

El  =  TEMP  ♦  SI  +  El  ♦  Cl 

0483 

TEMP  =  82 

0484 

82  =  82  ♦  C2  -  E2  ♦  S2 

0485 

E2  =  TEMP  ♦  S2  E2  ♦  C2 

0486 

TEMP  =  83 

0487 

83  =  83  ♦  C3  -  B3  ♦SS 

0488 

E3  =  TEMP  ♦  S3  B3  ♦  C3 

0489 

45 

TEMP  =  80  +  82 

0490 

82  =  80  -  82 

0491 

80  =  TEMP 

0492 

TEMP  =  81  +  83 

0493 

83  =  81  -83 

0494 

81  =  TEMP 

0495 

TEMP  =  BO  +  B2 

0496 

E£  =  BO  -  B2 

0497 

BO  =  TEMP 

0498 

TEMP  =  B1  +  B3 

0499 

B3  =  B1  -  B3 

0500 

B1  =  TEMP 

0501 

8Z 1=8 0+81 

0502 

8Z£=B0+B1 

0503 

8Z3=80-81 

0504 

8Z4=B0-B1 

0505 

8Z5=82-B3 

0506 

8Z6=B2+83 

0507 

8Z7=82+E3 

0508 

8ZS=B£-83 

0509 

997 

F0PM8T  (;iXf4f;E12.5»  1X>> 

0510 

8<:k:0>  =  CMPLX(8Z1»8Z£) 

0511 

8(:K1>  =  CMPLX  (8Z3!i  8Z4> 

0512 

8(:K2>  =  CMPLX  (8Z5»  8Z6> 

0513 

8(:K3>  =  CMPLX  (8Z7»8Z8> 

0514 

999 

FDRM8T  aX»'8  COEFFICIENTS^ 

0515 

50 

CONTINUE 

0516 

C 

WRITE  <37»998>  8  <K0>  » 8  (Kl>  » 8  CKE)  » 8  CKS) 

0517 

C 

WRITE  <37!. 995)  K;0»K1»KE»K:3 

0518 

995 

F0RM8T  <:1X!.4<;I4!.  1X>) 

0519 

998 

F0RM8T<:1X»4<;E12.5!.  1X.E12.5»  1X)> 

0520 

IF  <K  .LE.  1)  60  TO  55 
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’-I  OJ  ('■j  ^  U"'  'Xi  K  CO  «-■  CU  C'*  't  ll">  ■•£'  ^-  CO  '7'  O  •-<  (U  Ol  rf  Ip  -Xi  iDO  (7'  C»  «  (U  (*) 
fo  OJ  OJ  OJ  OJ  ru  OJ  OJ  OJ  00  oo  oo  co  oo  co  >00  00  00  00  ^  n-  tj-  rr  ir>  uo  uo  m 

uo  bo  u*i  u"'  b"i  uo  ir.i  bo  bO  10  IP  u"'  u"'  bo  bo  bO  bo  iT'  uo  uo  uo  bo  bo  bo  bo  bO  bo  bO  bO  bO  uo  bO  ip 


K  =  K  -  £ 

GD  TD  30 

55  KE  =  K:3  ISP 

IF  <KN  .LE.  KB>  GD  TD  70 
IF  CJ  .NE.  1>  GD  TD  60 
K.  =  3 
J  =  MK 
GD  TD  30 
60  J  =  J  -  1 
C£  =  Cl 

IF  <J  .NE.  £>  GD  TD  65 
Cl  =  Cl  ♦  CK  ♦  SI  ♦  SK 
SI  =  SI  ♦  CK  -  C£  ♦SK 
GD  TD  35 

65  Cl  =  <:C1-S1>  ♦  SQ 

SI  =  (:C£  ♦  SI)  ♦  SQ 
GD  TD  35 
70  CONTINUE 

C  WRITE  <37!.  398)  CFl  (KZ)  » KZ=1 ,  £56) 

IF  CM  .LE.  1)  GD  TD  9005 

MR  =  M  +  1 
JJ  =  1 
IWK<1)  =  1 
HD  75  I  =  £»MP 
IWKCI)  =  IWKd-D^E 
75  CONTINUE 

N4  =  IWK<MP-£) 

IF  <M  .GT.  £)  N8  =  IWK(MP-3) 

N£  =  iwk<;mp-i) 

LM  =  N£ 

NN  =  IWKCMP)  ♦  1 
MP  =  MP  -  4 
J  ==  2 
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0554 

80 

JK  =JJ  +  N2 

0555 

fiK2  =  8<:j) 

0556 

R<.j>  =  r<:jk:> 

0557 

R<:JK>  =  Rk2 

0558 

j  =  j-n 

0559 

IF  <;jj  .GT.  N4>  GD  TD  85 

056  0 

JJ  =  JJ  +  H4 

0561 

GO  TO  105 

056£ 

85 

JJ  =  JJ  -  N4 

0563 

IF  CJJ  .GT.  HS>  GD  TD  90 

0564 

JJ  =  JJ  +  MS 

0565 

GD  TD  105 

0566 

90 

JJ  =  JJ  -  N3 

0567 

K  =  MP 

0568 

95 

IF  tlWKCK)  .6E.  JJ>  GD  TD 

1  00 

0569 

JJ  =  JJ  -  IWKCKJ 

0570 

K  =  K  -  1 

0571 

GD  TD  95 

057£ 

1  0  0 

JJ  =  IWK<K>  +  JJ 

0573 

105 

IF  <JJ  .LE.  J>  GD  TO  no 

0574 

K  =  NN  -  J 

0575 

JK  *  hd  -  JJ' 

0576 

RK2  =  R<;j> 

0577 

R  <:  j>  =  R  <;  J  j> 

0578 

R  < JJ)  =  RK2 

0579 

RK2  =  R<K) 

0580 

RCK)  =  RCJK) 

0581 

RcJK)  =  RK2 

058£ 

no 

J  =  J  +  1 

0583 

IF  <J  .LE.  LM)  GD  TD  80 

0584 

9005 

2X=12 

0585 

CRLL  TYPECZX) 

0586 

C 

WRITE  <37.  998)  <.R  CKZ)  .  KZ=1 

.  256 

0587 

RETURN 

0588 

END 

0589 

C 

0590 

C 

0591 

SUBRDUT I NE  TYPE  <ZX) 

0592 

c 

WRITE <37. 20) ZX 

0593 

20 

FDPMRT</."  SUBROUTINE  ENTERED 

0594 

RETURN 

0595 

END 

:  \ 
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4.5 


COMPUTER  SOFTWARE  VERIFICATION 


To  assure  that  the  software  works  properly,  the  follow¬ 
ing  verification  procedure  has  been  established.  Before  running 
the  program  on  real  data,  it  is  recommended  that  the  program  be 
run  on  an  idealized,  sinusoidal  waveform  with  slight  distortion 
added . 

Two  cases  of  distortion  will  be  introduced;  the  first 
one  involves  amplitude  only,  and  the  second  one  involves  phase 
only.  The  distortions  are  presented  as  sinusoids  of  various 
amplitudes  but  always  small  enough  so  as  to  permit  the  use  of 
the  paired  echo  theory  concepts  to  be  valid.  Also,  the  period 
of  the  sinusoid  is  always  harmonically  related  to  the  bandwidth 
of  the  equalizer  network  and  as  such  is  consistent  with  the  con¬ 
cept  of  representing  the  distortion  as  a  Fourier  series. 

The  test  cases  of  distortion  are  separated  into  two 
groups,  one  of  which  assumes  only  amplitude  distortion  and  the 
other  of  which  assumes  only  phase  distortion,  but  of  varying 
levels.  The  FFT  presents  the  coefficients  for  the  taps  which 
in  effect  represent  the  amplitude  and  phase  of  the  echo  associ¬ 
ated  with  each  tap.  These  coefficients  are  then  compared  to  the 
calculated  values  based  on  the  paired  echo  theory. 

CASE  1:  A  sinusoidal  amplitude  distortion  of  various  peak  to 
peak  values. 

|G|  =  10  log^Q  n/i^+Q^  =  S.log^Q  (I^+Q^)  =  amplitude 
0  =  tan"^  ^  =  phase 
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To  simplify,  we  consider  the  effects  of  only  amplitude  distortions 
so  that  assuming  no  phase  distortion  and  letting  Q  =  0  we  get 

|G1  =  5  log^Q  =  10  log^Q  I 

0  =  tan"^  (^)=  0 

Since  a  sinusoidal  distortion  is  needed: 

sjn  =  10  log^Q  I 

where  AX  =  incremental  frequency 

N  =  number  of  points 

K  =  scale  factor  (for  K  =  1,  numerator  goes 
0  to  -2  or  2  dB 

and  10  log  I  represents  the  distortion  level  in  dB  ( peak-to-peak ) . 

In  order  to  fit  the  RADC  data  form  so  that  -128  <  I  <  128,  we  multiply 
the  given  values  of  I  by  128. 

TEST  1:  Amplitude  distortion,  0.5  dB  peak-to-peak  with  first 
echo  pair 

.0.5  =  *  10  log,Q  therefore  K  =  4 

For  215  points  in  2it  radians  or  360°,  there  AX  =  and  0  <N  <214* 

I  =  10  1A-X.NJ  ,-l] 

From  the  test  run  we  get  the  first  echo  peak  value  distortion  to  be 
-30.44  dB.  Analytically  (for  small  amp  distortion)  we  have: 

♦Note  that  the  number  of  sample  points  is  not  germane  to  the  analysis 
(as  long  as  it  is  a  large  number)  and  relates  only  to  the  expected 
number  of  samples  to  be  used.  We  later  change  this  to  256  points. 
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®  1 

20  IoQiq  ■'tude  ripple  in  dB  (peak  value) 


a,  1^1 

for  1  i  2  ^  ®c^’0  amplitude 

*0  ^  “o 


20  109io  (gTf)  “  echo  in  dB 
0 


20  log^Q  (1+  =  0.5  dB 


®1  10  [(f^)-i] 

j~  »  - - ^ -  =  0.0296 


echo  »  20  log^o  ^27“^  -30.566  dB 


TEST  2:  Distortion  is  introduced  in  the  second  echo  pair  with  an 
amplitude  ripple  of  0.5  dB 


I  »  10 


where  Ax  * 


2x360 

TTr 


and  0  <  N  <  215 


From  the  test  run  we  get  second  echo  peak  value  to  be  -30.53  dB  where 
the  calculated  value  is  the  same  as  before  -30.566  dB. 

TEST  3:  Distortion  is  introduced  in  the  first  and  second  echo  pairs 
with  an  amplitude  ripple  of  0.5  dB 


I  »  (10  [- 


sin  (Ax,N)-1 


-3  +  10  [- 


sin  (Ax,N)-l 


where  Ax^  -  |y|;  AX2  =■  and  0  <  N  <  214 
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Because  we  are  summing  two  sinusoidal  amplitude  distortions  of  0.5 
dB  each,  the  resultant  distortion  needs  to  be  divided  by  2  in  order 
for  it  to  be  0.5  dB . 

From  the  test  run  we  get  both  the  first  and  second  echo 
peak  to  be  -36.2  dB .  The  calculated  first  and  second  echo  peak  is: 


20  lo^lO  (1  "  4i-) 


10 


=  0.0148 


20  log  10  a 


=  36.58  dB 


Tests  4  through  6  are  the  same  as  Tests  1  through  3  respectively, 
except  the  peak  amplitude  distortion  is  set  at  1  dB  instead  of 
0.5  dB  (then  K  =  2).  Test  7  is  the  same  as  Tests  4  or  5  except 
the  distortion  is  introduced  into  the  eighth  echo  pair.  The 
results  are  given  in  Table  4-1. 

TEST  8:  Oistortion  is  introduced  in  the  first  and  eighth  echo  pairs 
with  an  amplitude  ripple  of  1  dB  (K  =  2) 


I  =  (10 


sin  (AX,N)-1  sin  (&X,N)-1 

•  A  _  t 


where  A'’  *  36^.  _  8x360 

wnere  a  ^  2i5>  ^^2  215 


and  0  <  N  <  21 4 


CASE  II;  A  sinusoidal  phase  distortion  of  various  peak-to-peak 
values.  As  defined  previously  the  gain  and  phase  are: 

|G|  »  10  log,Q  /  I^+Q^ 

9  =  tan”^  ^ 


Now  we  want  to  have  a  sinusoidal  phase  variation,  so  0  *  A  sin  (AXN)k. 
To  simplify,  let  us  set  |G|  =  0  dB  so  that: 
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TABLE  4-1.  TEST  RESULTS  OBTAINED  WITH  DEC-20  COMPUTER 


Test  1 
Test  2 
Test  3 
Test  3 
Test  4 
Test  5 
Test  6 
Test  6 
Test  7 
Test  8 
Test  8 


*Actua1  Distortions  outputted  by  the  algorithm  are  in  echo  pairs. 
The  value  shown  is  an  average. 
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1 /?  1 /? 

(1)  10  log^Q  [I^+Q^]  =  0  and  [I^+Q^]  =  1  or  I^+Q^  =  1 

(2)  j  =  tan  (A  sin(AXN)-k) 

Q  =  I  tan  (A  sin  (AXN)-k) 

Let  A  sin  (  XN)-k  =  0  so  that  Q  =  I  tan  0 
From  ( 1 )  we  ha ve : 

=  1  -Q^ 

From  (2)  tan^  9  =  (1-Q^)  tan^  0  =  tan^  0  -Q^  tan^  0 

P  P  9  P 

Q  +  Q  tan  0  =  tan^  0 

(1+tan^  0)  =  tan^  0 

Q  =  tan  0 

1+  tan  0 

o  1/2 

I  “  rl-  tan^  01  .  1 _ 

L  p  j  ^  ^  V/ 7 

1+tan^  0  (1+tan  0) 

0  =  tan”^  ^  =  phase 

To  simplify,  let  us  always  work  in  the  I-V  quadrant  so  that  Q  <  0 

- 2 — ’ 

and  I  >  0,  therefore  when  tan  0^0,  then  /  1  +  tan  0  >  0.  In  order 
to  fit  the  RADC  format  (2's  complement)  we  multiply  I  and  Q  by  128. 

TEST  1;  Sinusoidal  phase  distortion  of  6°  peak  value  is  intro¬ 
duced  in  the  first  echo  pair 

0  *  3  sin  (AXN)  -k  where  AX  »  0  <  N  <  214  and 

k  *  offset  *  42®  places  0  in  fourth  quadrant 
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and  I 


Q  =  tan(3  sin  (axN)-42) 
\l'[+  tan^(3  sin(AxN)-42) 


h 


From  computer  simulation  we  get  the  first  echo  to  be  -31.49  dB . 
Analytical ly  we  get 

27rx  b, 

’°9l0  (473^^  = 

where  b.|  =  ripple  in  degrees 

6,  .  6°,  20  log^,  (j^)-  -31.64  dB 

TEST  2:  Cosinusoidal  phase  distortion  of  6°  peak  value  is  intro¬ 
duced  in  the  first  echo 


0  =  3  cos  (AXN)  -42 


The  computer  simulation  is  basically  the  same  except  that  for  an 
odd  (sin)  distortion  function  the  phase  angle  for  the  two  echoes 
is  180°  apart. whereas  for  an  even  (cos)  distortion  function  the 
two  echoes  are  in  phase.  The  opposite  is  true  for  amplitude  dis¬ 
tortion.  That  is,  for  odd  function  distortions  the  phase  angle 
for  the  two  echoes  is  in  phase  whereas  for  even  distortion  function 
the  phase  angle  for  the  two  echoes  is  180°  apart. 

TEST  3:  Sinusoidal  phase  distortion  of  22.92°  peak  value  is 
introduced  in  the  first  echo 

0  =  11.46  sin  (AXN)  -42 
Calculated  result  is  20  log-jQ  ( 

simulated  result  is  -20.02  dB. 


The  results  for  the  phase  distortion  test  cases  are 
shown  in  Table  4-2. 
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TABLE  4-2.  TEST  RESULTS  OBTAINED  WITH  DEC-20  COMPUTER 


Peak  Phase 
Distortion 
( degrees ) 

Echo 

Simulated 

Di storti on 
Phase  (dB) 

Cal cul ated 
Di stortion 
Phase  (dB) 

Test 

1 

6 

1 

-31  .49 

-31  .64 

Test 

2 

6 

1 

-31  .49 

-31 .64 

Test 

3 

22.92 

1 

-20.02 

-20.0 

4.5  PROGRAM  TEST  ON  THE  HP  2100A 

Three  tests  were  run  on  the  algorithm  contained  in  the 
HP  2100A  computer  at  RADC.  The  first  test  is  a  single  echo  test 
involving  a  peak  sinusoidal  amplitude  distortion  of  0.5  dB.  The 
second  test  is  a  double  echo  test  involving  a  peak  sinusoidal  amp¬ 
litude  distortion  of  0.5  d3  with  one  sine  wave  at  twice  the  fre¬ 
quency  of  the  other.  The  third  test  is  a  single  echo  test  involv¬ 
ing  a  peak  phase  distortion  of  six  degrees. 

To  perform  the  first  test,  a  program  has  been  implemented 
(now  on  disk  cartridge  45)  called  "ATIO".  A  listing  of  this  program 
...is  gijvefl  in  Appendix  3. 

ATIO  outputs  M  sample  points  of  a  distorted  sine  wave 


given  by: 


I  =  10 


[sin  AXN]-1 
40 


where:  AX  =  and  N  =  0,1,2,...M 

In  our  case,  we  have  run  ATIO  with  256  sample  points  (M  =  255). 
As  mentioned,  the  phase  comoonent  of  this  test  is  0.  ATIO  pro¬ 
duces  both  the  single  echo  pair  and  two  pair  test  cases.  The 
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case  just  described  (one  echo  pair)  is  printed  out  in  a  data  file 
called  "ADI800";  this  file  contains  256  octal  mode  (06  integers, 
in  two  records,  128  numbers  each;  this  file  is  the  sequential 
output  of  the  "I"  or  amplitude  function  as  N  goes  from  1  to  256. 

To  run  the  double  echo-pair,  amplitude-distortion-only 
case,  we  executed  ATIO  again  but  such  that  ATIO  adds  to  the  first 
waveform,  a  second  one  given  by 

sintAX2N)-l 


where:  and  N  =  0,1,2, _ M  =  255 

The  superimposed  waveform  (this  is  what  is  printed  out  into  the 
data  file)  is 

■sup  =  “  "  ‘2>/2 

Again,  there  are  256  samples,  written  into  two  128-number  records 
in  06  format.  The  file  containing  is  called  "AD001".  This 

double  echo-pair  that  is  equivalent  to  Test  3  of  Table  4-1  which  shows 
the  theoretical  peak  distortion  to  be  -36.58  dB. 

For  the  third  case,  phase  distortion  only,  the  program 
"TEST  1"  was  executed.  The  waveform  produced  by  "TEST  1"  is  a  sine 
whose  phase  is  given  by: 

tan  (3  s  i n (  AX'!  N  ).42 ) 

^  ;  ’  X  =  11^,  N  =  0,1  ,2 . 255 

l+tan*"  (3  sin  (AXiN)-42) 

To  run  this  test  case,  we  executed  "TEST  1"  (also  on 
disk  45),  and  used  the  data  file  "AD002".  This  results  in  a  test 
case  of  6°  phase  distortion  equivalent  to  Test  1  in  Table  4-2. 
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Tables  4-3  through  4-5  give  printouts  of  the  data 
contained  in  the  output  files  for  these  three  test  cases,  "AD0180", 
"AD001",  and  "AD002"  respectively. 

We  can  now  compare  the  results  obtained  with  the  test 
cases  run  on  the  DEC-20  and  HP  2100A  computers  with  their  respec¬ 
tive  equalizer  algorithms.  Slight  differences  are  to  be  expected 
because  of  the  differences  in  precision  and  the  internal  mathematical 
algorithms.  Table  4-6  shows  the  results. 

The  comparison  shows  excellent  agreement  between  the 
two  equalizer  algorithm  outputs  (DEC-20  and  HP  2100A)  and  demon¬ 
strates  that  they  are  performing  the  way  they  are  supposed  to. 
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TABLE  4-3.  OUTPUT  TEST  1  (ONE  ECHO  PAIR) 
File  ADOOO 


Tap  No. 

Attenuator 

(dB) 

(Turns) 

Line  St 
(deg) 

retcher 

(Div) 

8A 

-.59754E+02 

9.924 

.77894E+03 

7.2 

7A 

-.57500E+02 

9.133 

.84826E+03 

7.8 

6A 

-.50899E+02 

9.931 

.72472E+03 

6.7 

5A 

-.54626E+02 

9.280 

.1 2426E+04 

1 1  .0 

4A 

-.76775E+02 

9.870 

.13916E+03 

1  .3 

3A 

-.57659E+02 

9.054 

.48888E+03 

4.5 

2A 

-.63654E+02 

8.213 

.82400E+03 

7.6 

1  A 

-.30575E+02 

2.894 

.94832E+03 

8.8 

0 

.95424E+01 

0.000 

.44309E+03 

4.1 

IB 

-.30102E+02 

2.571 

.40519E+02 

0.4 

2B 

-.58178E+02 

8.552 

.60589E+03 

5.6 

3B 

-.53561E+02 

9.354 

.47294E+03 

4.4 

4B 

-.57641E+02 

9.961 

.80104E+03 

7.4 

5B 

-.52656E+02 

9.388 

.15492E+03 

1  .4 

6B 

-.68077E+02 

9.716 

.23711E+03 

2.2 

7B 

-.55111E+02 

9.231 

.61405E+03 

5.7 

8B 

-.58767E+02 

9.938 

.91032E+03 

8.4 
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TABLE  4-4.  OUTPUT  TEST  2  (DOUBLE  ECHO  PAIR) 

File  ADOOl 


Tap  No. 

Attenuator 

(dB) 

(Turns) 

Line  Stretcher 
(deg) 

(Div) 

8A 

-.58976E+02 

9.646 

.78231E+03 

7.2 

7A 

-.58854E+02 

9.603 

.86500E+03 

8.0 

6A 

-.50607E+02 

9.816 

.85919E+03 

7.9 

5A 

-.55226E+02 

9.481 

.90270E+03 

8.4 

4A 

-.74326E+02 

9.505 

.32962E+03 

3.0 

3A 

- . 56551 E+02 

8.710 

.361 37E+03 

3.3 

2A 

- .36679E+02 

3.518 

.93857E+02 

8.7 

lA 

- . 36044E+02 

3.898 

.95012E+03 

8.8 

0 

.95424E+01 

0.000 

.44309E+03 

4.1 

IB 

-.35637E+02 

3.589 

.39872E+03 

3.7 

2B 

-. 35844E+02 

2.950 

.49132E+03 

4.5 

3B 

-.52810E+02 

9.094 

.60045E+03 

5.6 

4B 

-.57316E+02 

9.856 

.97058E+03 

9.0 

5B 

-.53165E+02 

9.563 

.13478E+03 

1  .2 

68 

-.66063E+02 

9.314 

.10264E+03 

0.9 

7B 

-.56190E+02 

9.654 

.59731E+03 

5.5 

SB 

-.58025E+02 

9.706 

.90695E+03 

8.4 
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TABLE  4-5.  OUTPUT  CASE  3  (PHASE  DISTORTION  ONLY) 

File  AD004 


Tap  No. 

Attenuator 

(dB) 

(Turns) 

Line 
( deg) 

Stretcher 

(Div) 

8A 

-.59834E+02 

9.953 

.79680E+03 

7.4 

7A 

-.59511E+02 

9.837 

.73877E+03 

6.8 

6A 

-.50643E+02 

9.831 

.96246E+03 

8.9 

5A 

-.56282E+02 

9.841 

.11 579E+04 

10.7 

4A 

-.72772E+02 

9.269 

.38686E+03 

3.6 

3A 

-.57695E+02 

9.066 

.39645E+03 

3.7 

2A 

-.66269E+02 

8.720 

.67455E+03 

6.2 

lA 

-.3n02E  +  02 

2.984 

.10384E+04 

9.6 

0 

.95424E+01 

0.000 

.44309E+03 

4.1 

IB 

-.31 104E+02 

2.74^ 

,17235E+03 

1  .6 

2B 

-.58022E+02 

8.500 

.44508E+03 

4.1 

3B 

-.53881E+02 

9.466 

.78023E+03 

7.2 

4B 

-.57201E+02 

9.818 

.86518E+03 

8.0 

5B 

-.53831E+02 

9.794 

.10693E+03 

1  .0 

6B 

-.68513E+02 

9.803 

.20575E+03 

1  .9 

7B 

-.56951E+02 

9.961 

.33624E+03 

3.1 

8B 

-.57922E+02 

9.674 

.66682E+03 

6.2 
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TABLE 

4-6. 

COMPARISON  OF  TEST  RESULTS  FOR  THE 
EQUALIZER  ALGORITHMS 

DEC-20  AND 

HP  2100A 

Test 

Type  of 

Di  stortion 

Echo 

Affected 

Simul ated 
Level  in 
HP  2100A 

Echo 

dB 

DEC-20 

Theoretical 
Level  in  dB 

Test 

^  , 

0.5  dB 

1  A 

-30.58 

-30.44 

-30.57 

One 

Echo 

Pair 

Peak 

Ampl i tude 

On!  y 

IB 

-30.10 

(AVG) 

-30.57 

Test 

2. 

lA 

-36.04 

36.36 

-36.58 

Ooub1 e 
Echo 

Pair 

0.5  dB 

IB 

2A 

-35.64 

-36.68 

(AVG) 

-36.58 

-36.58 

2B 

-35.84 

-36.58 

Test 

3, 

6°  Peak 

lA 

-31  .10 

-31 .49 

-31 .64 

Phase 

Distortion 

Only 

IB 

-31  .10 

(AVG) 

-31 .64 
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While  this  is  not  100%  conclusive,  it  does  provide 
a  high  degree  of  confidence  that  the  algorithm  is  working  properly 
The  only  real  test  is  the  use  of  the  algorithm  over  many  real 
situations  which  completely  exercise  it.  Some  comments  about 
what  to  expect  are  in  order. 

First  of  all,  the  Microwave  Transversal  Equalizer 
(MTE)  to  be  used  to  provide  corrections  to  the  distortion  has 
its  own  iherent  second  order  distortions  which  are  not  constant. 
That  is,  under  one  set  of  amplitude  and  phase  adjustments  for 
each  of  the  taps  there  is  a  given  "self-distortion"  which  is 
included  in  the  data  supplied  to  the  algorithm  to  analyze. 

When  the  amplitude  and  phase  settings  are  now  changed  to  cor¬ 
rect  for  the  total  distortion  reflected  in  the  data,  then  the 
new  settings  of  the  MTE  will  correct  much  of  the  distortion  but 
because  the  MTE  settings  are  different  from  the  original  settings 
there  will  be  a  new  "self-distortion"  introduced  by  the  MTE. 

This  means  that  the  correction  is  imperfect  and  the  process  must 
be  iterated  until  the  distortion  residue  is  less  than  the  allow¬ 
able  value  (or  until  the  "sel f-di stortion"  changes  are  greater 
than  the  corrected  indicated  by  the  algorithm). 

Another  source  of  distortion  lies  in  the  hardware 
used  to  provide  the  data  samples.  These  distortions  are: 

1.  the  phase  response  of  the  components 

2.  the  amplitude  response  of  the  components 

3.  noise  introduced  on  the  data  (Gaussian) 

4.  random  spikes  (transients  induced  by 
i nterference ) 
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5.  sinusoidal  variations  (coupled  interference 
from  motors,  oscillators,  etc.) 

6.  timing  errors  (clocks,  circuit  switching 
speeds,  etc.) 

7.  inaccuracies  in  data  samples  (improper 
scaling,  etc.) 

8.  digitizing  errors  (bit  errors  in  digital 
operations) 

These  errors  introduce  their  own  "self-distortion"  and  just  as  in 
the  MTE  limit  the  residual  distortion  that  can  be  obtained.  In 
addition,  some  of  the  errors  are  random  so  that  iterations  won't 
help. 

Finally,  there  is  the  algorithm  itself.  It  cannot  be 
considered  fully  debugged  until  all  combinations  of  operation  have 
been  exercised.  This  is  not  a  practical  thing  to  do  in  the  labor¬ 
atory  but  instead  requires  a  perigd  of  field  use.  In  addition,  the 
algorithm  has  limits  of  precision  which  will  ultimately  limit  the 
amount  of  distortion  that  can  be  corrected  even  if  the  hardware 
were  perfect,  although  this  is  the  area  that  is  least  likely  to 
present  a  problem. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Computer  algorithms  have  been  successfully  developed  to 
provide  open-loop  adaptive  control  of  the  MTE.  The  algorithms  are 
based  upon  the  application  of  FFT  techniques,  and  the  necessary 
corailary  software  programming  procedures  have  been  developed  and 
described . 

Verification  of  the  aforementioned  FFT  and  associated 
software  programming  has  also  been  accomplished  at  AIL  with  a 
dec-20  computer.  Such  verification  has  involved  the  oeneration 
of  output  adjustment  data  for  MTE  control  based  upon  computer 
analysis  of  artificially  simulated  time  sidelobe  input  distortion 
levels.  A  similar  procedure  will  be  undertaken  at  RADC  using 
real-time  I  and  Q  data  acquired  at  the  A.  Froelich  High  Power 
Tube  Facility.  Such  distortion  levels  will  be  processed  at  RADC 
using  translation  of  format  programs  between  the  AIL  DEC-20  and 
the  RADC  HP  2100A  computers.  Preliminary  data  suggests  that  the 
post-delivery  program  translation  will  be  successful  after  the 
normal  debugging  procedures  have  been  completed. 

AIL  recommendations  for  future  MTE  related  efforts 
include  the  following: 

0  Continuation  and/or  extension  of  contracted 
efforts  to  include  AIL  post-delivery  support 
for  RADC. 
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- - - - 

®  Investigation  and  definition  of  suitable  closed- 
loop  techniques  and  procedures  for  adaptive  oper¬ 
ation  of  MTE. 

•  Development  of  suitable  solid-state  time  delay 
devices  to  replace  the  present  manually  adjust¬ 
able  line  stretchers. 

•  Modification  of  the  MTE  to  expand  capability  by 
providing  32  tap  operation. 

•  Development  of  suitable  electronic  interface  to 
implement  fully  adaptive  closed-loop  operation 
of  the  MTE. 
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APPENDIX*  1 

SAMPLE  COMPUTER  RUN  WITH 
INPUT  DATA 
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APPENDIX  1:  SAMPLE  RUN 


This  appendix  includes  a  sample  run  of  the  program  with  the  input  data 

entitled  AD000. 


:  RU«  HIL 
1 

INITIhL  sett I mgs  RT  step  0  SHOULD  BE: 


HTTENUHTOP 


stpptchep 


1  RslO.OCiO 

2  Rsl 0. 000 

3  R=1  0.000 

4  R*1  0.000 

5  0.000 

6  R=1 0.000 

7  R=1  0.000 
S  R= I  0.000 

M  *  0. 000 

1  B*1  0.000 

2  B*l 0.000 

3  B=1 0. 000 

4  B*1  0.000 

5  B=1  0.000 

6  E=1 0.000 

7  B=l  0.000 

8  B»10. 000 


1 

R* 

7.3 

R* 

G.  2 

3 

R= 

V  •  t 

4 

R= 

1.5 

•• 

R= 

9.  1 

Initial 

R= 

7.5 

Attenuator 

»■’ 

R= 

7.  0 

Settings 

8 

R= 

6.  1 

M  = 

4.  1 

1 

E= 

1.2 

ci 

B= 

3.7 

3 

B= 

5.2 

4 

B= 

7.2 

5 

B= 

.5 

6 

B= 

1.4 

7 

B= 

3.2 

3 

B= 

6. 2 

Initial 

Line 

Strecher 
Sett  inns 
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APPENDIX  1  SAMPLE  RUN  (continued 
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APPENDIX  1  SAMPLE  RUN  (continued) 


INTERMEDIRTE  SETTINGS  RT  STEP  1  RRE: 

RTTENURTOR  LINE  STRETCHER 


1 

R= 

1 . 376 

1 

R= 

O  ^  C; 

a 

R* 

6.  303 

3 

R= 

1*  •  C* 

3 

R= 

3.513 

3 

R= 

li.  5* 

4 

R= 

ci.  436 

4 

H= 

1 . 6 

5 

R= 

3 . 663 

5 

R= 

3. 6 

6 

R= 

3.705 

6 

R= 

3.  3 

R= 

3.  333 

7 

R= 

3.  3 

3 

R= 

3.  4  34 

3 

R= 

6*1 

M  = 

0.  000 

M  = 

4.  1 

1 

B= 

1.573 

1 

.4 

c 

B- 

6 .  035 

E= 

5-5 

3 

B= 

3.  665 

3 

B= 

6 .  0 

4 

B* 

3.  363 

4 

B= 

:  c  1 

5 

E= 

3.  707 

5 

B= 

1 .  0 

6 

B* 

3.  334 

6 

B= 

■Z*  m  Z* 

l' 

B* 

3.  343 

"7 

B= 

3.  3 

s 

B= 

3.503 

3 

B= 

3  •  c! 

RFTER  0  ITERRTIONS  DISTORTION  CRN  NOTtBE  CONPENSRTED 
flIL  :  STOP 


Note:  XD  *  20  in  $AT10  when  generating  input  for  this  run. 

In  subsequent  runs,  XD  «  40 
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APPENDIX  2 

PROGRAM  TO  GENERATE 
TEST  DATA 
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APPENDIX  2.  PROGRAM  TO  GENERATE  TEST  DATA 
This  routine  generates  ideal  sinusoidal  distortion  as  the  data  file  for  the  oroaram 


0  0  0 1 

FTH4. 

L 

0  ij  Oil' 

PPGbPRM  RT 1 0 

0  0  0  5 

DIMENZIOR  IEUF3<35b>  «  IEUF3:  •:35b''  .NEUF':35b.' 

Ij  Ij  ij4 

D I  MEN  ;  I  ON  1  DC  E  ■  3?3  •  .  NRM  <  3  •  •  I  ECF  0  >:.35bJ  •  I  3  1 3E  (.3 

0  'j  ij5 

D I  MEN  ;  I  ON  LEUF  ..35b  >  .  I EUF  1  <  35b  ■ 

ijijije' 

D I  MEN  7  I  □r4  I  DC  E3  •;.373  •  •  NRM 1  3  :• 

1  j  0  0  r 

IT', 'PE  =3 

0  0  ijy 

i:i3£=3 

0  0  O'? 

I  3IZE  •.3>  =133 

Ij  0 1  0 

DR TR  rtRM  3HRD  •  3H 0  0  •  3H 0  ' 

0  0 1  1 

DRTR  NRMl  "3HhD« 3H00* 3H1  ■ 

0  0 1  c 

CRLL  CF'ERl  IDCE.  lERP.NRM,  ICIZE*  ITYPEJ 

Ij  0  I  7: 

CRLL  CF'ERT  <  I  DC  E !  .  I  EF'P  •  NRM  1  .  I  '3 1 ZE »  I  T YF'E> 

0  0  1  4 

ri  =  35b 

Ij  Ij  t 

:;ri=4  0. 

0 1  j  1  4 

D:<1  =3bO.  •  FLORT  vN-l  j 

IJ  'J  1  7 

DX3=L,/:i>3. 

I.I  0  1  : 

DO  .'-,  0  .1=1 'N 

Ij  0 1  ■? 

:;=D:a»Fi.0RT*:  j-i' 

j  ijc  Ij 

'C=Ii.''.3  ♦FlOhT  'ij-l  ■ 

’  j  Ij  i  1 

1  =  '  1  0.  ■  IN  1415.''130.  '•  -1  >  .'XD>  > 

0  0?i 

;'3=  <  1  0.  ;  IN  <Y>3.  1415 '130.  >  -1,>  'XD  '  > 

Ij  iji  :i 

LEUF  '  .|•  =  133♦/:l 

'jOi'4 

XM=,  .'l-t-X3-  '3. 

ijiji'5 

MEi  IF  ,  J:'  =  133^XM 

OOib 

I F  ■  MEUF  c  j  ;•  .  i5E .  1 33  •  MEUF  j:.  a  137 

0  OcT' 

MEUF  <  J  ■•  =35b^MElJF  •;  J> 

OOciS 

IF  <LEiJF  <  J  '  .EQ.  1  337' LEUF  < J  '  =  1 37 

OOi"? 

LEijF  (  J.'  =35b'*’LEUF  (.  J> 

0  0  ?  0 

30 

CONTINUE 

0  Ij  5 1 

DO  30  1=1 '133 

0  0  Sc' 

lEUFO  jl  ■  =LBUF  a  '■ 

0  Ij  3  5 

IEUF3  <  I  •  =MEUF  I  I  :• 

00  34 

30 

CONTINUE 

00  3  J 

DO  4  0  k  =  1  .  1  33' 

0  0  3b 

I  EUF  1  ( k  :•  =LEUF  <  1 33-*-K  ’ 

0  0  3  ? 

I  EUF  3 • k ■ =MEUF  < 1 33+k  > 

0  033 

4  0 

CONTINUE 

00  3'? 

CRLL  MR  I TF  a  DC  E • I ERR • I EUF  0 • 

0040 

C  RLL  MR  I  TF  <  1  DCE  •  I  ERR  -  I  EUF  l  :• 

004  1 

ChLL  MRITF  < IDCEl . IERR« lEUFE' 

0043 

CRLL  i-iR  I  TF  •:  I  DC  E 1 .  I  ERR »  I  EUF  3  • 

004  5 

CRLL  CLOIEaDCE- lERR' 

0  044 

CRLL  CL0:E  » IDCEl. I  ERR' 

0045 

'TOP 

004b 

END 
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This  program  generates  a  sinusoid  with  phase  distortion  only. 


LIST.STESTl: :45 

■STESTl  T=00004  IS  DM  Cfi'00045  USING  00005  ELKS  R=0034 
FlVLSLl  FTN4 


nOOii 
0  0  0  5 
0  0  04 
0005 
0  0  06 


PPDGRRN  TESTl 

INTEGER  I  BCE  <£?£>  »NRM(:S-f>  <•  IEUP0<£5G>  » lEUFl  tESG;-  .  LEUF  (£56;- 
INTEGER  ISIZE<:£> 

ITYPE=£ 

N=£5G 


0007 
00  OS 
00  09  C 
'  010 
0  0 1 1 
0  0 1  £ 

0  0 1  3 
0  0 1 4 
0  Cl  1  5 
0  Cl  1 6 
0  0 1  7 
Cl  0 1  Si 
0  0 1 9 
0  Oc'  0 
00£1 
CiCi££ 

00  £3 
00  £4 
00  £5 

00£b 
0  0  £7 
00£S 
00  £9 
0  0  3  0 
0  0  3 1 
0  Ci3£ 

0  Cl  3  5  1  0 

00  34 
8035 
00  3  8 

0  0  3  7  3  Cl 

0038 
0039 
004  0 

0041 


IS:IZE=£ 

ISIZEC£>=1£S 

CRLL  CPERTariCE!.  lERR.NRM,  ISIZE.  ITYPE) 
DRTR  NRM.'  £HRD«  £H00»  £H4 
CRLL  CRERTCiriCE.  IERR»NRM,  IS12E»  ITYPE> 
DY 1 =36  0 .  FLORT  CN- 1 > 

DO  10  I=1.N 
:<:i=Dxi  ♦FLORT  a-i;' 

X  =  3 .  ♦  S I N  •:  X 1  ♦  3 .  1 4 1 5  1 8  0 .  >  -4£ . 
X£=S1N(X^3.  1415.'180.) 

X3  =  CDS  <X^3. 1415/180. ) 

T=X£/X3 

T£=T^T 

X4=S0RT  (:i-*-T£  :' 

IFCT.GT.O)  X4=-X4 
Q=T/X4 
I  Q=  1  £3 .  ♦1.1! 

XI=SClRT<:i.-Q^Q> 

II=1£8^XI 

IFCII.GE. 1£7>  II=1£7 

IF  (II.LE.-lESl-  II='1£8 

IF  tIQ.GE.l£7:'  1 0=1  £7 

IF  CIO. LE. -1£8>  I0=-l£8 

IFCII.LT.  O:-  II  =  II+£56 

IF  (IQ.LT.O-  IG!=IQ+£56 

LEUF<:n=£56^II>IQ 

CONTINUE 

DO  30  J=1j1£8 

IB  UFO  Cj:-  =LBUF 

I EUF 1 ( J>  =LEUF  < J+ 1 £8> 

CONTINUE 

CRLL  WRITFCIDCB. IERR» IBUFOJ 
CRLL  URITFUDCB.  lERRj  IBUF1> 

CRLL  CLOSE  UDCE*  lERR) 

STOP 
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Input  for  Case  1:  Octal  Format 


:LIST?ftD00i;i 

ft  DO  00  T=000  0c:  IS  DM  CROOOO 


RECs:  00001 

074000  074400  074400  074400 
075000  075000  075000  075000 
075400  075400  075400  076000 
076000  076000  076400  076400 
076400  076400  077000  077000 
077000  077000  077000  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077400  077400  077400  077400 
077000  077000  077000  077000 
076400  076400  076400  076400 
076000  076000  076000  076000 
075400  075400  075400  075000 
075000  074400  074400  074400 

REC«  0000£ 

074000  074000  074000  074000 
073400  073400  073400  073000 
073000  072400  072400  072400 
072000  072000  072000  072000 
071400  071400  071400  071400 
071400  071400  071000  071000 
071000  071 0 0 0  0 710 0 0  0 7100 0 
071000  0710  0  0  0  7 1 0  0  0  0  7 1 0  0  0 
071000  071000  07 1 0 0  0  07 1 0 0 0 
071000  071000  0710  0  0  0  7 1 0  0  0 
071000  071 0  0  0  0  7 1 0  0  0  071 0  0  0 
071400  071400  071400  071400 
072000  072000  072000  072000 
072400  072400  072400  072400 
073000  073000  073000  073000 
073400  073400  073400  074000 


USING  00002  BLKS  R=01£8 


074400  074400  074400  075000'< 
075000  075400  075400  075400^ 
076000  076000  076000  076 000* 
076400  076400  076400  0764 00* 
077000  077000  077000  0770 00^ 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077400  077400* 
077400  Ci77400  077400  077400* 
077400  077400  077400  077400* 
077400  077400  077000  077000* 
077000  077  000  077000  07‘64'j0* 
076-400  076.400  076000  076' 000* 
076000  075400  075400  075400* 
075000  075000  075000  075000* 
074400  074400  074400  074000* 


074000  073400  073400  073400* 
073000  073000  073000  073000* 
072400  072400  072400  072400* 
072000  072000  072000  072000* 
071400  071400  071400  071400* 
0 7 1  0  0  0  0 7 1  0  0  0  0  7 1  0  0  0  0  7 1  0  0  0  * 
0  7 1  0  0  0  0  7 1  0  0  0  0  7 1  0  0  0  0  7 1  0  0  0  * 
0 7 1 0  0  0  0  7 1 0  0  0  071 0  0  0  0  7 1 0  0  0  * 
071000  071000  071000  071000* 
0  7 1 0  0  0  0  7 1 0  0  0  071 0  0  0  071 0  0  0 ♦ 
071000  071400  071400  071400* 
071400  071400  071400  072000* 
072000  072000  072000  072400* 
072400  072400  073000  073000* 
073000  073400  073400  073400* 
074000  074000  074000  074000* 
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Input  for  Case  2;  Octal  Format 

LIST.flnOCil 

FtriOi;il  T=OOOOc:  is  oh  CROOOOc:  USIHG  0000c:  BLKS  R=01£8 
PEC"  00001 


074000 

074400 

074400 

0744  0  0 

075000 

075400 

075400 

075400 

076000 

0  7 6  0  0  0 

0764  0  0 

076400 

077000 

077  0  0  0 

077  000 

077  0  0  0 

077000 

077  0  0  0 

077400 

077400 

077400 

0774  0  0 

077400 

077000 

077  0  0  0 

077000 

077000 

077  0  0 1'l 

076400 

07ii^4ijij 

076400 

076400 

076000 

076000 

076000 

075400 

075  0 1'l  0 

075000 

C  ~  5  0  0  0 

075000 

074400 

074400 

074400 

074400 

074000 

074000 

074  000 

074000 

073400 

073400 

07:3400 

073400 

073400 

073400 

073400 

0734  0  0 

07:3400 

073400 

073400 

073400 

074  0  0  0 

074000 

074000 

074  000 

PEC”  0000£ 

074000 

074000 

074400 

074400 

074400 

074400 

074400 

074400 

075000 

075  0  0  0 

075000 

075000 

075  0  0  0 

075000 

075000 

075000 

074400 

0744  ijO 

074400 

0744  ij(i 

074400 

074400 

074  000 

074000 

074000 

073400 

0734  0  0 

073400 

073000 

073:000 

073000 

075  0  0  0 

07£400 

07£400 

07£400 

07£000 

07£000 

071400 

071400 

0714-00 

071400 

071400 

071400 

071400 

071 000 

071 000 

071 000 

071 000 

071400 

071400 

071400 

071400 

071400 

07£000 

07£000 

07£000 

07£400 

07£400 

07£400 

07£400 

073400 

073400 

073400 

073400 

0744  0  0 

075000 

075000 

075 000^ 

075400 

076000 

076000 

076000-* 

076400 

0764  0  0 

076400 

076400^ 

077000 

077000 

077000 

077 000^ 

077400 

077400 

077400 

0774 00^ 

077000 

077000 

077000 

077000* 

077000 

077000 

0764  00 

0764 00* 

076400 

076000 

076000 

076000* 

075400 

075400 

ij754  0  0 

075400* 

075000 

075000 

074400 

0744 '! 

0744  0  0 

074000 

074000 

074  000* 

074000 

074000 

074000 

074 OoO* 

073400 

07:54  00 

073400 

073400* 

073400 

0734  0  0 

073400 

075400* 

073400 

074000 

074000 

074000* 

074000 

074000 

074000 

074000* 

074400 

074400 

074400 

0744  0  0* 

0744  0  0 

074400 

074400 

075000* 

075000 

075000 

075000 

075000* 

075000 

075000 

075000 

075000* 

074400 

074400 

074400 

074400* 

074000 

074000 

074000 

074000* 

073400 

073400 

073400 

073000* 

073: 0  0  0 

07£400 

07£400 

07£4  0  0* 

0  7£  0  0  0 

07£000 

07£000 

07£000* 

071400 

071400 

071400 

071400* 

071 000 

071 000 

071 000 

071000* 

071 000 

0710  0  0 

071 000 

071400* 

071400 

071400 

071400 

071400* 

07£000 

07£000 

07£000 

07£400* 

073000 

075000 

073-:000 

073000* 

074000 

074000 

074000 

074000* 
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Input  for  Case  3:  Octal  Format 


L I ST • hD  0  04 

FIDO 04  T=0000c:  IS  ON  CROOOOc  USING  0000c'  ELKS  R=01£S 


REC”  00001 

057653 

057653 

057653 

057653 

057654 

060 £5 4 

06  0£54 

060£54 

060c 55 

06  0£55 

o»;s*=; 

OA  riA®i5 

06  0656 

060656 

060656 

060656 

06 1 £56 

061 £57 

061 £57 

061 £57 

061 £57 

061 £57 

061 £57 

061 £57 

061660 

061660 

06 1 66  0 

06 1660 

06 1 66  0 

061660 

06 166  0 

06 1 66  0 

061660 

06 1660 

06 166  0 

06 1 66  0 

06 1 66  0 

06 1 66  0 

06 1660 

06 166  0 

06 1 66  0 

061657 

06 1 657 

06 1657 

061 £57 

06 1 £57 

061 £5 7 

061 £57 

06 1 £56 

06 1  £56' 

061 £56 

060656 

060656 

060656 

06  0655 

06  0655 

060£55 

060 £55 

06  0254 

060254 

OF' 

057654 

057654 

057654 

057653 

REC«  0000£ 

057653 

057£53 

057£53 

057252 

057£5£ 

057£5£ 

056652 

05665 1 

056651 

flSAA®!! 

Ci5665 1 

05665 1 

056£50 

056£50 

056250 

056250 

055647 

055647 

055647 

055647 

-QF^  coii£--e«> 

055647 

055647 

055647 

055647 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055£46 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055246 

055646 

055647 

055647 

055647 

055647 

055647 

055647 

056250 

05625  0 

056250 

056250 

056250 

056250 

056251 

05665 1 

056651 

05665 1 

056651 

057252 

057252 

057252 

057252 

057653 

0576.53 

057654 

057654« 

060254 

060254 

06  0255 

060255< 

060655 

060655 

06  0655 

060656* 

06.0656 

061256 

06 1 £56 

061256« 

061 £57 

061257 

061257 

061£57» 

06'  1 257 

061657 

061657 

06 1 66  0^ 

06 1 6.6-  0 

06- 1 6.6. 0 

06 166  0 

ij6  1 66. 0^ 

06 166  0 

06 1 66  0 

0 6 1 6  6  0 

06. 1 66  0^ 

06 1 6-6  0 

06 1 66  0 

U  t'  1 1 1*  LI 

06 1 66  0^ 

06 1 660 

06 1 66  0 

06 1 66  0 

i  ft-  0^ 

06- 1  £57 

061257 

061257 

IJ6  1257*- 

061257 

061257 

061257 

1 

060656 

060656 

06  06.56 

06  0656^ 

06  0655 

060655 

06  0255 

060£55^ 

060254 

060254 

06  0254 

060254^ 

057653 

057653 

057653 

057653^ 

057252 

057252 

057252 

057£52» 

05665 1 

056651 

056651 

05665 !♦ 

05625  0 

056250 

l.l  5 5  L* 

•55625  •5* 

056250 

Ci56c'5Ci 

056250 

0.55647» 

055647 

055647 

055647 

055647^ 

055647 

055246 

055246 

055£46^ 

055246 

055246 

055246 

i555246^ 

•555246 

055246 

055246 

055£46« 

055246 

055246 

055246 

•555£46^ 

055246 

055246 

055246 

055246* 

055647 

•555647 

•555647 

•555647  ♦ 

•555647 

•555647 

055647 

•5556.47  ♦ 

056250 

•55625  ^5 

•55625^5 

•55625  •■•♦ 

•556651 

•556651 

•556651 

05665 !♦ 

056652 

•557252 

I557252 

057252« 

057252 

•557253 

•557653 

•557653^ 
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1 


Input  for  Case  3: 
Decimal  Format 


24491 

24491 

24491 

24491 

24491 

24491 

24492 

24492 

24492 

24745 

24748 

24748 

2474?. 

24743 

24749 

24749 

24749 

2474? 

25!X<5 

2'C<!W 

25iVj5 

25005 

25.;0.6 

25006 

‘7'^r‘O^, 

25006 

250r6 

252i2 

2526: 

252S2 

25262  ■ 

'A  ^ 

'  ■*  ,  A  ■i 

■*S7A.? 

-c  ".2  “ 

25263 

25519 

25519 

:5'7'. 

^520 

25520 

25520 

25520 

25520 

25520 

25.520 

1552>’i 

25"2'.'- 

2:520 

‘-S2.". 

25520 

25520 

25520 

25520 

25520 

23520 

25520 

25.320 

25520 

2>5‘20 

2-52' ) 

25520 

25'20 

23520 

25520 

25520 

25520 

25519 

f  3 

2551? 

25263 

252S3 

25263 

25263 

25:63 

25263 

25263 

'5C.*,  5 

252c3 

25262 

^  ~  *.'/ 

«‘'V. 

25.:«j6 

25005 

25005 

25005 

25005 

24749 

24749 

24749 

24749 

2474? 

24748 

2474A 

2“" ‘♦'I* 

24'4? 

2474? 

24492 

24492 

24492 

24491 

2449' 

24491 

24491 

24491 

24491 

242'35 

24235 

24234 

24234 

24234 

24234 

:'42'94 

24224 

2423“ 

23975. 

23977 

79,g77 

23977 

23977 

23977 

22.977 

23977 

23977 

:3.720 

'.7720 

■:3720 

23720 

23720 

23720 

23720 

23720 

23720 

2.3726 

23720 

2:'hS3 

22-463 

23463 

22463 

23463 

23463 

224-3 

204-;:'? 

23463 

23463 

2'-.463 

234.63 

23.463 

23206 

2:2i>S 

23206 

23206 

23206 

23206 

23206 

23206 

23206 

to  •'*>  w 

23206 

2.3206 

23206 

23206 

2320S 

23206 

23206 

23206 

23206 

2.3206 

2320'S 

237,'. 6 

25206 

to  :i.vo 

23206 

232s:« 

23206 

23206 

■23206 

23206 

2:.106 

23206 

23206 

23206 

23462 

23463 

23463 

23463 

23^63 

23463 

23463 

23469 

23463 

r<463 

234.93 

234  63 

■23463 

23463 

23720 

23720 

23720 

23720 

23 '/"20 

'■3720 

23720 

23720 

23720 

23.721 

23977 

23977 

23977 

23977 

23977 

23977 

23977 

23977 

2397.8 

24234 

24234 

24234 

24234 

24234 

24234 

24234 

24234 

24235 

24491 

24491  STOP 

I 


I 
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MISSION 

of 

Rome  Air  Development  Center 

MVC  plan6  and  txtaUu  AeAeMch,  dtveZopmtnt,  tzAt  and 
AzZccteA  acqiuAition  paogfiam  In  Aixppofvt  o(,  Command,  ContAot 
CommunicationA  and  Inteltigzncz  IC^i)  aaUvitiu.  TzcknicM 
and  zngX.nzzAA.ng  AuppoAt  uiithin  oazoa  tzchnizaJt  zompztzncz 
AM  pKooldzd  to  ESP  PAogAom  O^^XczA  (POa)  and  othzA  ESP 
ztmzntA.  Thz  pAtnctpal  tzzhnical  nuAAton  oA-Zoa  oaz 
communccat<.onA,  ztzcXAomagnztic  giUdancz  and  zontAoZ,  Aua- 
vzitZanzz  o^  gAound  and  azAoApacz  obizctA,  tntzZCtgzncz  data 
cottzctcon  and  handting,  tn^oAmation  AyAtzm  tzzhnotogy, 
tonoAphzAic  pAopagation,  Aotid  Atatz  AztznczA,  mtcAouuvz 
phyAtcA  and  zlzctAonic  Aztiabttcty,  matntainabtUty  and 
compattbtCtty. 


END 


DATE 
FILMED 


DTIC 


